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ABSTRACT 



The objective of the STAR project is to test and space 
qualify a new continuous cycle cryogenic refrigeration 
system for cooling of sensors and electronics which is based 
upon the newly discovered thermoacoustic heat pumping 
effect. The new refrigerator has no sliding seals, a cycle 
frequency of about 300 Hz, and uses acoustic resonance to 
enhance the overall power density and efficiency. This 
thesis is concerned specifically with the design and testing 
of the electrodynamic transducer which is responsible for 
the electro-acoustic power conversion. A computer model of 
the driver/resonator system is presented along with the 
techniques for measurement of the electrical and mechanical 
parameters used as input for the model. 

A final driver design (including dimensional drawings) 
utilizing a modified JBL 2450J neodymium-iron-boron 
compression driver and associated leak-tight electrical 
feed-throughs, microphone, accelerometer, pressure gage, 
pressure housing, and resonator interface is provided. 
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INTRODUCTION 



I . 



A. BACKGROUND 
1 . History 

Thermoacoustics can generally be described as the 
study of the interaction between heat and sound. Scientific 
interest in this interaction is not new. Lord Rayleigh (1878 
and 1945) discussed various qualitative examples of the 
production of sound oy heat. In one of these examples he 
describes investigations by Sondhauss (1850) of an effect 
observed by glassblowers : heating a bulb of glass at the 

end of a hollow rod occasionally led to the generation of 
sound. Faraday (1818) demonstrated that this effect occurred 
with gases other than air. Sondhauss proved that the vibra- 
tion of the glass itself did not generate the sound, but he 
offered no explanation as to what did. In his description of 
Sonhauss ’ work, Rayleigh stressed the importance of the 
phase difference between temperature and particle motion. 
Rijke (1859) built an open tube with a wire mesh inside. 
When the mesh was heated, the tube produced a sound of 
impressive intensity. The functioning of this oscillator is 
slightly different since it requires steady gas flow for its 
operation. A similar sound production effect was observed by 
Taconis , et al. (1949) in hollow tubes immersed in liquid 

helium ( ^ ) . The Taconis oscillations were addressed 
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quantitatively by Yazaki , et al . (1980). The work described 
above deals primarily with prime movers— devices that con- 
vert a temperature gradient to sound energy. We call this 
the classical era of thermoacoustics. 

Interest in thermoacoustics was renewed when the 
idea of the use of acoustical devices as refrigerators 
appeared. Gifford and Longsworth (1966) described the pump- 
ing of heat along a surface caused by a periodic change in 
the pressure of the adjoining gas. Such a change can be 
produced by the oscillations of a sound wave. In their 
experiments Merkli and Thomann (1974) explored the heating 
and cooling effects on the wall of a gas-filled resonant 
tube. They found that heat was transported from a region 
near the velocity antinode (or maximum value) of the sound's 
standing wave to regions near the adjacent pressure anti- 
nodes. Thermoacoustic prime movers have also received recent 
attention. Kempton (1976) investigated the excess noise of 
aeroengines above that predicted by theory. He determined 
that the sound was produced by unsteady heat transfer. Each 
of these latter three groups of exprimenters used some 
theory for comparison, but it was mostly qualitative. None 
provided the complete quantitative theory that would explain 
their experimental observations. 

The theoretical breakthrough in thermoacoustics 
field was made by Nicklaus Rott (1969, 1974, 1975, and 
1980). His theory combined basic principles from physics, 



thermodynamics, and acoustics to quantitatively describe the 
effects found in both types of thermoacoustic devices: prime 
movers and heat pumps. Prime movers use a temperature gra- 
dient to create sound, as in the effects discussed by Sond- 
hauss, Rayleigh, Taconis , and Kempton . Heat pumps, or refri- 
gerators, use the oscillating pressure of a sound wave to 
produce a temperature difference, as described by Gifford 
and Longsworth and Merkli and Thomann. Rott described the 
effect found by Merkli and Thomann as thermoacoustic stream- 
ing. It is this effect that makes thermoacoustic refrigera- 
tion possible. 

Inspired by Ceperley's (1979) traveling wave Ster- 
ling Cycle heat engine and Rott's quantitative theory, the 
team of Wheatley, Hofler, Swift, Migliori , and Garrett 
(1982, 1983a, 1983b, 1985, and 1986) developed a series of 
thermoacoustic experiments at Los Alamos National Laboratory 
in New Mexico. They investigated the basic thermoacoustic 
effects in both prime movers and refrigerators and compared 
their experimental results to Rott’s theory. 

2 . Thermodynamics 

We'll digress from history here to explain the ther- 
modynamic distinction between prime movers and refrigera- 
tors. This discussion follows Sears and Salinger (1975). A 
prime mover receives heat at a high temperature, T M (the hot 
reservoir), does work on its surroundings, and rejects heat 
at a lower temperature, T c (the cold reservoir) , as shown in 



Figure I-l.a. The first law of thermodynamics tells us that 
the heat flow (Q M ) from the hot reservoir must equal the 
work done on the surroundings (W) plus the heat flow (Q c ) to 
the cold reservoir: 

Q M = W + Q c . (1-1) 

The second law of thermodynamics states that in every pro- 
cess the entropy of the universe must either remain constant 
or increase. The entropy change for the isothermal heat 
transport to/from the reservoirs is defined as the heat flow 
divded by the reservoir temperature. For a prime mover the 
second law can be written: 

^ m i iSg or Qh/Th a Go/To. ( 1 — 2 ) 

The thermal efficiency of the prime mover is defined as the 

ratio of the work output to the heat input: 

n = W/Qm = (Oh - Qc ) /Oh = 1 - Qc/Qh . (1-3) 

Comparing equations 1-2 and 1-3 gives us a limit on the 

efficiency based on the temperatures of the reservoirs: 

h — 1 — T c /Tm. (1-4) 

This limit is known as the Carnot efficiency. 

In a refrigerator the process is basically reversed 
(see Figure I-l.b). Heat (Q c ) is removed from the low tem- 
perature reservoir, work (W) is received from the surroun- 
dings, and heat (Q M ) is delivered to the higher temperature 
reservoir. The expression for the first law is the same as 
for the prime mover, Equation 1-1 . For the refrigerator it 
means that the heat flow into the hot reservoir must equal 
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(b) 



Figure 1-1. The Thermodynamics of (a) Prime Movers 
and (b) Refrigerators. 
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the heat flow from the cold reservoir plus the work done to 
the system. The second law for this case is: 

ASm £ ASc or Qh/Th £ Qc/Tc . (1-5) 

The efficiency parameter for the refrigerator is called the 

coefficient of performance (COP) . It is defined as the ratio 

of the heat flow from the cold reservoir to the work input 

to the refrigerator, or: 

COP = Qo/W = Oc/(Q h -Qc) £ Tc/(T„-Tc). (1-6) 

The limit is known as the Carnot coefficient of performance. 

The work in this thesis deals solely with j. efrigera- 
tors. Even though prime movers are mentioned from time to 
time, the primary focus of the rest of our discussion will 
be refrigerators. 

3 . Overall Efficiency 

According to Hofler (1986), the early thermoacoustic 
work at Los Alamos focused on experimental refrigerating 
engines. The performance of these engines fell short of 
expectations, leading to simple experiments on basic thermo- 
acoustic effects and a proof-of-principle experiment on 
thermoacoustic refrigeration. Hofler then applied the Rott 
theory to the experimental systems and solved the resulting 
equations numerically. For his doctoral dissertation from 
the University of California, San Diego, Hofler designed and 
constructed a completely functional thermoacoustic refriger- 
ator. He also made accurate measurements of its thermody- 
namic efficiency, and used this efficiency to make compari- 
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sons to the Rett theory. After receiving his doctorate, 
Hofler came to the Naval Postgraduate School (NPS) as a 
post-doctoral fellow, and brought his prototype refrigerator 
with him. 

The purpose of this thesis, in conjunction with the 
work of several other students, is to modify Hofler’s re- 
frigerator design in order to improve its overall efficien- 
cy and make it suitable for space cryocooler applications. 
The basic design of the thermoacoustic refrigerator is shown 
schematically in Figure 1-2. The driver (A), which produces 
the sound waves, is coupled via bolts to the resonator (D) 
via a reducer cone (B) and bellows (C) . Inside the resonator 
is a stack of plastic plates (E) and their heat exchangers 
(F) which allow the heat to be removed from the hot end and 
absorbed by the cold end. It is the interaction between the 
sound waves and the plastic plates that produces a tempera- 
ture difference across the plates and/or pumps heat from the 
cold heat exchanger to the hot heat exchanger. A brief 
description of efficiency would be useful here. Figure 1-3 
shows a diagram of the chain of energy transformations which 
occurs in a thermoacoustic refrigerator. There are three 
conjugate variable pairs (V and I; P and U; and AT and 6 0 ), 
and two transformations ( electrical- to-acoustical and acous- 
tical-to-thermal ) . 

We start with electrical input power to the driver, 
which can be calculated (using r.m.s. values) with: 
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Figure 1-2. Basic Thermoacoustic Refrigerator Showing 

(A) Driver, (B) Reducer Cone, (C) Bellows, 
(D) Resonator, (E) Stack, and (F) Heat 
Exchangers . 
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Figure 1-3. Chain of Energy Transformations in a 
Thermoacoustic Refrigerator. 
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p, 



1 •<= 



i v cos e, 



(1-7) 



where I is the current into the electroacoustic driver, V is 
the voltage across the driver terminals, $ is the phase 
angle between I and V, and cos £ is known as the "power 
factor." An electrodynamic driver converts the electrical 
power to acoustic power (the first transformation) . Analo- 
gous with the previous definition of the electrical input 

power, the acoustic power is given by: 

W = P U cos 0 , (1-8) 

where P is the acoustic pressure, U is the volumetric velo- 
city which equals the product of particle velocity and 
resonator cross-sectional area, and $ is the phase angle 
between them. The efficiency of this electroacoustic trans- 
formation is given by : 

= ft/P EI _ EC . (1-9) 

In the second transformation, the refrigerator con- 
verts the acoustic power to a temperature gradient (vT) and 
heat flow (6). As discussed previously, the efficiency for 
this transition is given as the coefficient of performance: 

COP = (5c/«. (1-10) 

The overall thermoacoustic efficiency is therefore the pro- 
duct of the efficiencies for the two transformations: 

n-r« = IIea * COP = Qc/Peuec . (I-ll) 

Improvements in the overall refrigerator efficiency 

can thus be divided into two distinct, though interrelated, 

efforts. The first is the subject of this thesis: modifica- 
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tions to the driver to optimize the electroacoustic 
efficiency. The second effort involves the optimization of 
the coefficient of performance and is addressed by LT 
Michael P. Susalla, USN , in his master’s thesis. 

B. SPACE THERMOACOUSTIC REFRIGERATOR (STAR) 

1 . Motivation 

In addition to improving overall efficiency, our 
experimental refrigerator will be designed and built with 
the intention of flying it on the Space Shuttle as part of 
the National Aeronautics and Space Administration’s (NASA) 
Get Away Special (GAS) program. As of 3 February 1988, there 
is a Memorandum of Agreement between the Naval Postgraduate 
School (NPS) and the Air Force which provides funding for 
the space flight and which assigns the NASA payload, number 
G-337 to this project. The reason for testing the thermo- 
acoustic refrigerator in space is that the immediate appli- 
cations of the STAR are expected to be the cooling of elec- 
tronics, high T c 

superconductors, and infrared detectors in 

space . 

There are currently two primary methods for cooling 
satellite systems in space: (1) evaporation of expendable 

cryogens (liquid helium, nitrogen, ammonia-methane, and 
solid hydrogen, etc.), and (2) closed cycle refrigerators 
(Stirling cycle, Vuilleumier cycle, etc.), which involve 
large reciprocating piston masses operating at low 
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(1984)] The 



frequency. [Walker (1983) and Smith, et al . 
disadvantages of these current cooling methods are their 
short lifetimes (expendable cryogens) and high vibration 
levels and low reliability (closed cycle refrigerators). The 
advantages of the thermoacoustic refrigerator over these and 
other systems lies in its simplicity and reliability. The 
STAR has no sliding seals, efficient performance, low vibra- 
tion levels, and (hopefully) a longer life span. 

The thermoacoustic refrigerator needs to be surroun- 
ded by a vacuum and insulation material to minimize heat 
loss due to thermal conduction and radiation. The vacuum of 
space will provide the best insulation to determine the 
refrigerator's ultimate efficiency. The absence of gravity 
removes the possibility of thermal convection as an addi- 
tional nuisance heat transport mechanism. 

2. Get Away Special (GAS) Program 

NASA’s GAS program allows small, self-contained 
payloads to fly on the Space Shuttle in GAS canisters for 
relatively low cost ($10,000) (Get Away Special Team, 1984). 
The GAS can is five cubic feet in volume and can house a 
payload of up to 200 pounds. Each payload must contain its 
own electrical power, control, data acquisition and storage 
facilities, etc. The Shuttle's astronauts will turn a switch 
on or off at designated times during the flight, but are 
otherwise not involved with the experiment. 
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3. 



Shared Subsystems from NASA G-313 



Another group of NPS students and faculty (Boyd, 
et al . , 1987) have taken advantage of the GAS program to 
measure the resonant acoustic modes of the shuttle payload 
bay and the ambient acoustic environment produced as a 
result of main engine and booster operation during launch. 
This experiment is titled "The Space Shuttle Cargo Bay 
Vibroacoustics Experiment" and is designated by NASA as 
payload G-313. Several subsystems that were developed for 
NASA G-313 will be used by NASA G-337 (STAR) . A schematic of 
the STAR in its GAS can is shown in Figure 1-4. One of the 
systems borrowed from NASA G-313 is the electronics system 
used to run the experiment and record the data . The recorder 
system consists of the INTEL model BPK 5V75 magnetic bubble 
memory module and an NSC 800 microprocessor-based control- 
ler. Two other NPS students, LT Charles B. Cameron, USN and 
CPT Ronald Byrnes, USA, will be designing the analog elec- 
tronics and software to integrate these systems into the 
STAR experiment for their master’s thesis. A block diagram 
of the electronics is shown in Figure 1-5. 

Another borrowed system is the power supply, which 
consists of Gates brand lead acid battery cells (five 
ampere-hour, two volts each). These gelled electi'olyte bat- 
teries are ideal for the STAR due to their high power den- 
sity, low cost, and the absence of outgassing during dis- 
charge cycles. NASA G-313 used a one layer battery of 68 
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cells providing 680 watt-hours of energy and weighing about 
80 pounds (including the cells' support structure). We will 
be using two battery layers with as many cells as we can, up 
to 136. This will give us a minimum of 680 watt-hours and a 
maximum of 1,360 watt-hours of available electrical energy, 
and a battery weight between 80 and 160 pounds . The number 
of battery cells we can use will depend on the total weight 
of the driver-resonator assembly and its auxiliary equipment 
(vacuum can, gas reservoir, etc.) and electronics. 

The use of the GAS can imposes certain restrictions . 
Since we are using batteries to supply the power, the refri- 
gerator has to be energy efficient. Also, the GAS can setup 
requires the STAR to be compact and lightweight. Figure 1-6 
shows a photograph of Hofler’s prototype refrigerator. This 
setup is approximately six feet high. In comparison, the 
maximum payload height for the GAS can is 28.25 inches, or 
less than 2.5 feet. These considerations played a major role 
in the choice of equipment for and the design of the STAR. 

C . SCOPE 

Chapter II discusses basic acoustical and thermodynamic 
theory as applied to STAR. The results of this theory deter- 
mine the driver requirements. The theory of the selected 
electrodynamic driver is given next. Following the driver 
theory is a presentation and discussion of a computer 
program that models the coupled system, which consists of 
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Figure 1-6. Photograph of Hofler’s (1986) 
Prototype Refrigerator. 
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the electrodynamic driver and the resonator containing the 
gas mixture and stack . 

Chapter III describes the technique for measuring the 
electrodynamic driver parameters. These parameters are used 
as input for the driver design, the computer model, and the 
efficiency calculations. This measurement technique was 
developed using an Altec 290-16K electrodynamic driver and 
was applied to the final driver. 

Chapter IV describes the measurement of the STAR driver 
parameters. The driver chosen for STAR is a custom unit, 
built by Harmon- JBL ! <TM> , which uses neodimium-iron-boron 
(NIB) magnets. The STAR driver is a minor modification of a 
new line of compression drivers using NIB magnets and has 
characteristics similar to the JBL 2450 J. 

Chapter V discusses the construction of the STAR driver, 
including the housing and accessories (reducer, stiffener, 
bellows, accelerometer. microphone, electrical feed- 
throughs, etc.) necessary to integrate it with the resonator 
and the GAS cannister system. Dimensional shop drawings for 
the parts are included in Appendix C. 

Chapter VI gives conclusions and recommendations for 
further development. 
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II . 



THEORY 



A. THERMOACOUSTIC THEORY 

Thermoacoustic theory has been developed in detail by 
Rott (1969, 1974, 1975 and 1980) and adapted to the thermo- 
acoustic refrigerator by Wheatley et al. (1982, 1983a, 
1983b, 1985 and 1986), Wheatley and Cox (1985), Hofler 
(1986) and Swift (1989). 

The space thermoacoustic refrigerator (STAR) basically 
consists of an acoustic driver producing sound waves in a 
resonant tube (see Figure 1-2). This resonator is filled 
with a mixture of helium-xenon gas (12.5% xenon) pressurised 
to ten atmospheres, and contains a stack of plastic plates. 
This Chapter will present a qualitative model for the ther- 
moacoustic heat pumping process followed by a quantitative 
development for heat and work flow at a plate. The Chapter 
concludes with the efficiency of a stack of plates that are 
very much shorter than one quarter of a wavelength. 

1 . Qualitative Picture 

Consider a parcel of gas that moves back and forth 
along one of the plates at the acoustic frequency (see 
Figure II-l). As it moves, the gas parcel will experience 
changes in temperature and volume. Part of the temperature 
changes come from the adiabatic compression and expansion of 
the gas by the sound pressure, and part as a consequence of 
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the local temperature of the plate itself. A temperature 
gradient may develop along the plate as a result of the 
operation of the refrigerator. The temperature and volume 
changes can be described by six separate steps (the changes 
in volume are indicated by the change in size of the square 
parcel in the figure) . 

Assume the plate is at a uniform temperature, T. The 
parcel starts at position 1 (X = 0) with a temperature of T. 
The sound wave moves the parcel to the right to position 2. 
The parcel has undergone adiabatic compression and its tem- 
perature is now T ++ . Since the temperature of the gas is 
higher than that of the plate, heat (Q) will flow from the 
gas parcel to the plate. The parcel's volume is decreased 
and its temperature lowered to T + . This heat flow also 
causes the plate’s temperature to increase at the position 
X + . Parcel position 3 is actually in the same location as 
position 2, but it is displaced vertically in the figure for 
clarity. The parcel now moves to position 4 and expands 
adiabatically to a new temperature of T- A repea t of this 
process puts the parcel in position 5 with a temperature of 
T . Now the temperature of the gas is lower than that of 
the plate and heat (Q’ ) flows from the plate to the gas 
parcel, expanding the parcel and raising its temperature to 
T~~. This flow causes the plate's temperature to decrease at 
position X- . The parcel moves to the right again under 
adiabatic compression and we are back where we started in 
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position 1 with temperature T (adapted from Wheatley, et 
al., 1985). 

This gas parcel cycle is repeated at the resonator's 
operating frequency. For the STAR this frequency is between 
approximately 250 and 600 cycles per second, depending on 
the gas mixture (specifically the speed of sound in the gas 
mixture) and the length and shape of the resonator. 

Notice that the temperature of the gas parcel at 
position X = 0 is different depending on which direction the 
gas parcel is moving. It xs this phase shift in temperature 
relative to motion that produces the thermoacoustic effect, 
as we will show in the calculations that follow. 

It is also important to note that in order for heat 
to flow between the gas parcel and the plate, the parcel 
must be vertically located within about a thermal penetra- 
tion depth of the plate. The thermal penetration depth (Sk) 
is the distance that heat can diffuse through the fluid 
during a time l/u>, where tu is the acoustic angular 
frequency . 

Wheatley et al . (1986) describe thermoacoustic en- 
gines as consisting of long trains of these gas parcels, all 
about a thermal penetration depth from the plate. The par- 
cels draw heat from the plate at one extreme of their oscil- 
latory motion and deposit heat at the other extreme. Adja- 
cent heat flows cancel except at the ends. The net result is 
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that an amount of heat Q is passed from one end of the plate 
to the other. 

2 . Single Plate 

Swift (1989) illustrated the basic principles of the 
thermoacoustic engine by examining a simple example of a 
single plate in a gas-filled tube such as we described 
above. In this example the acoustic and thermodynamic ef- 
fects are nearly distinct. In the absence of the plastic 
plate, the plane standing wave sustained in the tube has 
well-known acoustic properties. It’s pressure and x-direc- 
tion velocity (along the length of the tube) are given as 
(using peak values): 

P x = p« sin(x/A) ( I I — 1 > 

and u x = i(p«/p m a) cos(x/A) (II-2) 

where is the pressure amplitude at the pressure antinode, 
the reduced wavelength A = A/2rt = a/a>, A is the wavelength, 
a is the speed of sound in the gas, w is the angular fre- 
quency, and p m is the mean density of the gas. The subscript 
1 indicates the first order time oscillation, and the sub- 
script m indicates the mean value. The total pressure can be 
given as 

P = Pm + Pie 1 ^. ( 1 1 - 3 ) 

A similar expression represents the total velocity. Follow- 
ing Swift, we will assume that the first order in the acous- 
tic amplitude is sufficient for all acoustic and thermo- 
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dynamic variables (pressure, velocity, temperature, density, 
and entropy) . 

The sound wave is adiabatic and has an oscillatory 
temperature given by: 

Tt = (T m e/p m c„) Pi, (II -4 ) 

where $ = - (dp/dT),=./p m is the isobaric thermal expansion 

coefficient and c P is the isobaric (constant pressure) heat 

capacity per unit mass. Note that T x and Pi are in phase. 
For ideal gases, thermodynamics tells us that 

Tm 3 / P mCp — (u — l)Tfii/)iPfn, ( I I “5 ) 

where t - the ratio of isobaric to isochoric specific heats 
(5/3 for monatomic gases, smaller, but greater than one. for 
other gases). Combining Equations II-4 and II-5 gives us : 

T x /T m = [(#-1)/*] pi/p m . ( 1 1 -6 ) 

The introduction of the plastic plate into the 
standing wave modifies the original, unperturbed temperature 
oscillations. This modification is due to the heat flow 
between the gas and the plate. as described previously. The 
temperature is modified in both magnitude and phase for gas 
about a thermal penetration depth away from the plate. 
According to Swift, this results in two important effects: 
1) a time-avei age heat flow near the surface of the plate, 
along the direction of acoustic vibration, and 2) the 
generation or absorption of real acoustic power near the 
surface of the plate. In refrigerators the effect is an 
absorption of acoustic power. In prime movers the acoustic 
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power is generated by the temperature gradient present in 
the plate. 

There are several assumptions we will make to 
simplify calculations, again following Swift. We assume that 
the plate is short enough compared to a reduced wavelength 

(Ax « A) and far enough from both velocity and pressure 
nodes that p x and u x can be considered uniform over the 
entire plate. We assume the gas has zero viscosity, so that 
ui does not depend on y (we had already made this assumption 
by looking at u x in the x-direction only) . We assume that 
the plate has a large enough heat capacity per unit area 
that its temperature does not change appreciably at the 
acoustic frequency. We assume that the plate has a mean 
temperature gradient in the x-direction ?T m . Finally, we 
neglect the plate’s and gas’ thermal conductivity in the x- 
direction . 

Applying these assumptions, we see that the mean 
fluid temperature (T™(x) ) is the same as that of the plate. 
Swift calculated the oscillating fluid temperature using 
the general equation of heat transfer. He kept only first- 
order terms, neglected thermal conduction along x, and 
applied the boundary condition T x (0)=0 imposed by the plate. 
The resulting equation is: 

T i (y ) = [ ( Tm8/ PnCp )Pi-( vT m /sB) Ui ] [ 1-exp [~(l+i)y/S] ] . ( 1 1 -7) 
To interpret this equation, we look at it in the limit that 
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the gas is far enough from the plate (y»SK) to make negli- 
gible thermal contact with the plate. This gives: 

Tx (y)-K (T m 0/p mC ^)p l -(vT m /a J )u i ] (II-8) 

The first term in Equation II-8 is the temperature 

oscillation due to the adiabatic compressions and expansions 

of the gas (see Equation II-4). The second term comes from 

the mean temperature gradient in the gas . As the gas 

oscillates along x with displacement amplitude ui/oj, the 
temperature at a given point in space oscillates by an 

amount vT m Ui/.u even if the temperature 01 a given piece of 
fluid remains constant. The actual temperature oscillations 
are a linear superposition of these two effects. 

The y dependent part of equation I 1-7 is complex. It 
approaches 1 for y » Sk, as given previously. It approaches 
zero for y « Sk, where the plate imposes the condition 
Ti =0. Most importantly, for y =s Sk , its magnitude is still 
of the order of 1, but it has a substantial imaginary part. 
This phase shift in the oscillating temperature of the 
standing wave at y 2* Sk, due to the thermal presence of the 
plate, is an important result because it leads directly to 
the time-averaged heat flow in the x-direction. It is the 
same phase shift we emphasized in the description of the gas 
parcel motion. 

a. Heat Flow 

Swift argues that since we are neglecting ordi- 
nary thermal conductivity in the x direction, the only way 
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heat can be transported along x is by the hydrodynamic 
transport of entropy, carried by the oscillatory velocity 

U i : 

9= = T m Pm s x u i — 1/2 p mCpp Im [ T i ] Ui . (II— 9) 

The line above the quantity SxUx indicates the time-average 

of the first-order entropy and velocity product. The heat 
flow is a second order quantity, signified by the subscript 
2 . 

The total heat flow along the plate, in the x- 
direction, is found by integrating q 2 over the y-z plane: 

= IT 0 /^ 9 = dy (11-10) 

where IT is the perimeter of the plate in the y-z plane. 

Substituting for T t and performing the integration gives: 

= -1/4 ITSK(Tm0)PxUx(r-l) (11-11) 

where ITSk is the thermodynamically active area in a plane 
perpendicular to the longitudinal acoustic motion, T m 0 is 
the heat parameter of the gas, and F is defined as the ratio 
of actual temperature gradient to the critial temperature- 

gradient ( vT/vT«=r- 1 1 ) . The critical temperature gradient 
occurs when the temperature change along the plate just 
matches the temperature change due to the adiabatic 
compression of the gas, and no heat flows between the gas 
and the plate. It is the boundary between the refrigerator 
and prime mover functions of the thermoacoustic engines. 

Equation 11-11 shows that when f 1 , heat flows 
up the temperature gradient from cold to hot and work 
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(acoustic power) is absorbed, as for a refrigerator. When 
r = 1 there is no heat flow. When T > 1, heat flows down the 

temperature gradient from hot to cold and acoustic power is 
produced, as for a prime mover. 

Note that the total heat flow is proportional to 
the area ITSk, and to (=1 for ideal gases). It is also 

proportional to the product p l u i , and so equals zero if the 
plate is at either a pressure node or a velocity node of the 
standing wave. The maximum value of the product occurs 
hai ,-ay between the nodes. Finally, the heat flow is 
proportional to the temperature gradient factor T-l. For 
vT m > 7 T crlt ., F-l >0 and the heat flow is toward the pres- 
sure node, while for vT m < 7T C n«, T-l < 0 and the heat flow 
is away from the pressure node. If suitable heat exchangers 
at temperatures T M and T c are installed at the ends of the 

plate (with T M - To = vT m Ax) , this heat flow carries heat 
from one exchanger to the other. 

The heat flow is small under ordinary 

circumstances. However, in closed resonators it is possible 
to achieve sound amplitudes many orders of magnitude higher 
than those of ordinary conversation . Since <S = is 

proportional to p^ 2 , and since in practical acoustic engines 
the entire cross section of the standing wave is filled with 

plates (spaced roughly 4 Sk apart) , very high heat flows 
and/or a large T M - Ti= may be achieved. 
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b. Work Flow 

The work flow (i.e. acoustic power) is given by 
the work per cycle times the rate at which that work occurs 
(the acoustic frequency f). From thermodynamics, the average 
acoustic power produced per unit volume is: 

= - (a i/fm) ipipi — - 1/2 uiB pi Im [ Ti ] . ( 1 1 — 1 — ) 

The gas about a thermal penetration depth away 
from the plate "breathes," because of thermal expansion and 
contraction, with the right time phase with respect to 
oscillating pressure to do (or absorb) net work. This is 
exactly the same gas that we have seen is responsible for 
the heat flow. Gas elsewhere is ineffective in doing (or 
absorbing) work. The density oscillations for v ^ g K anc j *■ 

y » Sk are in phase with the pressure oscillations, and 
hence do (or absorb) no net work. 

The total acoustic power produced is found by 
integrating &= over all space, as with the heat flow: 

W 2 = 1/4 ITSk U-l) ( Pi 2 /? m a ) ( T - 1 ) Ax/?-., (11-13) 

where all the terms have been defined previously. 

The acoustic power is proportional to the volume 
ITSkAx of fluid that is about a thermal penetration depth 
from the plate. It is proportional to p t 2 , and so is quad- 
ratic in the acoustic amplitude (as was the heat flow) and 
vanishes at pressure nodes. Finally, W 2 j_ s proportional to 

(T-l), the same temperature gradient factor as appeared in 
(Equation 1 1- 11). When vT m = (T-l) = 0, and there 
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are no temperature oscillations in the fluid other than 
those due to adiabatic processes , and no acoustic power is 

absorbed or generated. For V Tm > vT =r - lt , r -1 > 0 and acous- 
tic power is produced near the plate. Whether this power 
increases the amplitude of the standing wave, is radiated 
away to infinity, is simply absorbed, or flows through an 
acoustic- to-electric transducer to generate electric power 
depends on details of the resonator, not on the plate itself 

or on the standing wave near the plate. For vT m < 7T c n»,, T- 
1 < 0 and acoustic power is absorbed near the plate. For a 

tube without plates which has a diameter less than the 
wavelength, at constant temperature ( vT = 0), this work flow 
is responsible for the ordinary thermal attenuation of a 
sound wave (Kinsler, et al., 1932). 
c. Efficiency 

We combine Equations 11-11 and 11-13 to get the 
efficiency of the plate with no viscous or longitudinal 
conduction losses: 



W =/G= 


= U-1)/(T„,6) (uj Ax px ) / ( p m a 2 ui ) . 


(11-14) 


o sin 


x/A, and 


Px = p m a uo cos x/A, 




n = 


U-1)/(T, 


n.0 ) ( Ax/A tan x/A) . 


(II-15) 


For x 


« A, n 


= U-l)/(T m 0) (Ax/x), 


(11-16) 




and 


n = Hcaii-not/r 


(II-17) 



This efficiency approaches the Carnot efficiency 
as the power output and heat transfer rates approach zero. 



30 



We can make a similar calculation for the 
refrigerator mode of operation, where the relevant 
efficiency is the coefficient of performance, COP = 

We find that 

COP = r COP^„^„ ot . (11-13) 

3 . Short Stack 

After developing the basic principles of the thermo- 
acoustic engine using a simplified example, Swift developed 
a more realistic model by including viscosity, longitudinal 
thermal conductivity, finite (instead of infinite) plate 
heat capacity, and many plates. He still made the assump- 
tions that the temperature spanned is much less than the 
absolute temperature, that the length of the plates is much 
shorter than the reduced wavelength, and that the spacing 
between the plates is greater than a few penetration depths . 

Using the equation of motion for the gas, the boun- 
dary condition that the velocity at the gas-plate interface 
is zero, the continuity equation of the gas, the heat flows 
in the gas and plates, and basic thermodynamic properties, 
Swift derived a wave equation for pressure P t (x) in terms of 
dT m /dx, material properties and geometry. Swift also derived 
an equation for enthalpy flow along x in terms of pi(x), 
T m (x), geometry, and material properties (including Prandtl 
number of the gas cr = c p ij/k, where u is the viscosity, and k 
is the thermal conductivity of the gas). 
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In the boundary layer approximation, for a short 
(Ax « A) stack of plates, and neglecting viscosity, the heat 
flux and acoustic power are given by expressions very simi- 
lar to the simple results of the single plate example. 
Practical engines can be expected to have all available 
cross-sectional area filled with plates spaced about 43 k 
apart. The only effect of longitudinal (i.e., along x) 
thermal conductivity is to add to the heat flow in the x 
direction by simple conduction. For gases this effect is 
normally negligible. 

The inclusion of viscosity adds considerable 
algebraic complication, and a little conceptual difficulty 
as now Ui is a function of y. 

Without viscosity, vT = nt. is the temperature 
gradient for which the temperature oscillation at a given x- 
location is zero. But in the presence of viscosity, T x 
depends on x as well as y, so no well-defined vT crlt . exists. 

To the lowest order in 'fa (where cr is the fluid’s 
Prandtl number), the hydrodynamic heat flow is increased by 

1/(1 -’/o’). This factor arises because the mean velocity <u x > 
underestimates the velocity with which entropy is convected 
along the stack. As we saw in the single plate calculations, 
the convective entropy transport occurs mostly at a distance 
ok from the plate. The velocity there is higher than the 
mean velocity by l/(l-ycf). 
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The lowest order effect of viscosity on the acoustic 
power is the power dissipated by viscous shear in the fluid 
within the boundary layer, a well-known fluid mechanical 
result . 

For accurate calculations, the full general equation 
must be used; but for rough estimates, the results presented 
here are good approximations, and are much simpler to 
compute with. In addition, the results presented here are 
fairly easy to interpret physically. The expressions for 
arbitrary viscosity are complicated, but those to lowest 
order in viscosity are simple. It is easy to see the single 
plate expressions for heat flow and acoustic power about a 
thermal penetration depth from the plates, modified as 
longitudinal thermal conductance adds to the heat flow, as 
viscous shear dissipates acoustic power, and as non-ideal 
plate properties modify the thermal boundary condition at 
the fluid-plate interface. 

4 . Design Considerations 

One of the important considerations in the design of 
a thermoacoustic refrigerator is the location of the stack 
of plates with respect to the pressure and velocity nodes 
and antinodes of the acoustic standing wave in the resonant 
tube. Experiments performed by Wheatley, et al . (1986) 
showed that the heat flow as a function of position fit a 
simple sine curve whose spacial period is half the wave- 
length of the acoustic standing wave. By noting how the sign 
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of the temperature difference varied with respect to the 
plate's location in the sound wave, they saw that heat 
always flows in the direction of the closest pressure anti- 
node. This effect is expected because a parcel of gas moving 
in the direction of a pressure antinode is compressionally 
warmed and will transfer heat to the plate. A parcel moving 
toward a pressure node is cooled by expansion and will draw 
heat from the plate. At both the pressure nodes (velocity 
antinodes) and antinodes (velocity nodes) heat flow drops to 
zero. Thus the acoustic heat flow depends on both the acous- 
tic pressure and the fluid velocity. 

As a plate or stack of plates is moved away from the 
pressure antinode, the temperature gradient developed 
becomes smaller. At a quarter of a wavelength, no gradient 
forms (the critical temperature gradient equals zero). This 
positioning effect is important in the design of a refriger- 
ator, because, together with the length of the plates, it 
places an upper limit on the maximum temperature drop 
possible across the stack. 

B. DRIVER REQUIREMENTS 

In order to drive the thermoacoustic refrigerator we 
needed a device, a transducer, that would convert the elec- 
trical energy of a battery pack to acoustic energy, as 
discussed in Subsection I. A. 3. We wanted as high a pressure 
(or force) and volume velocity (or motion) as possible to 
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get the most acoustic power provided to the system. There 
were several types of electroacoustic transducers to choose 
from: electro-static, electret, electrodynamic, magneto- 
strictive, electro-strictive , and piezoelectric. Only one of 
these, the electrodynamic transducer, can provide our system 
with both the large forces and velocities necessary. This 
transducer will be described in the next section. The 
remainder of this section describes the rest of the trans- 
ducers and explains why each was not appropriate for use in 
the thermoacoustic refrigerator. 



The 


electrostatic 


transducer 


consists of 


a 


pair 


of 


charged 


electrodes. 


or capacitor 


plates , one 


of 


which 


is 


held stationary while 


the other , 


the diaphragm, 


moves 


in 


response 


to electrical 


excitation . 


The electret 


is 


a kind 


of 



electrostatic transducer that has a diaphragm of polarized 
plastic, and therefore does not need an external voltage 
supply to the diaphragm (Kinsler, et al , 1932, p. 350). If 
an alternating voltage is applied across the plates, the 
diaphragm moves, thereby radiating an acoustic wave. 
Unfortunately there is very little force (or pressure) and 
displacement associated with this motion. This type of 
transducer may be useful as a microphone but is not an 
appropriate choice to drive the thermoacoustic refrigerator. 

In a magnetostrictive transducer a change in the 
magnetic polarization in the material causes an elastic 
strain (Wilson, 1985, p. 3.). Although the forces are more 
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than adequate the strain produces too little motion to be 
useful to us . 

In electrostrictive and piezoelectric transducers an 
externally applied electric field causes a change in the 
dielectric polarization of the material, which in turn 
causes an elastic strain (Wilson, 1985, p. 3). Once again 
there is sufficient force but the displacement due to the 
strain is too small to be useful. 

Because these techniques have sufficient force, it would 
be possible to overcome the displacement limitation using 
mechanical leverage or other forms of mechanical impedance 
transf ormation , but this was not considered due to the 
associated increase in device complexity. 

C. ELECTRODYNAMIC THEORY 

1 . Physics 

The electrodynamic driver is a moving coil 
transducer which converts electrical energy to kinetic ener- 
gy (motion) . It consists of a diaphragm attached to a cylin- 
drical coil of wire (voice coil) that is suspended in a 
fixed magnetic field (see Figure II-2) . If an alternating 
current is supplied to the coil, the interaction of the 
current and the magnetic field will induce a force on the 
coil so that the diaphragm moves. The magnitude of the force 
is F = B1I, where B is the magnitude of the magnetic induc- 
tion, 1 is the length of the coil and I is the current. The 
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Figure II-2. Schematic of an Electrodynamic Transducer 
Showing (A) Magnet, (B) Back Plate, 

(C) Pole Piece, (D) Acoustic Seal, 

(E) Diaphragm, (F) Voice Coil, 

(G) Surround (Suspension) , and 

(H) Electrical Lead Connection. 
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direction of the force. given by the cross product of the 
current vector and the magnetic field (from Ampere's law F = 
II X B) , is everywhere orthogonal to the current and the 
radial magnetic field. The motion of the diaphragm produces 
sound waves with a frequency equal to that of the alterna- 
ting current. Another important property of the electrodyna- 
mic transducer is that the motion of the voice coil induces 
an electromotive force (emf) which equals B1 times the 
velocity of the voice coil (emf = Blv) . This "back emf” 
increases the electrical impedance of the driver and can be 
used to monitor the driver motion as discussed in Section 
III . B. 

The driver is attached to the resonator via a 
reducer and bellows. The reducer is a cone-shaped piece of 
aluminum mechanically connecting the driver voice coil which 



produces 


the 


force to 


the bellows. which 


in 


turn is 


in 


contact 


with 


the gas 


in the resonator. 


The 


reducer 


is 



designed to be lightweight, rigid, and strong for a direct 
transfer of motion and energy from the moving coil to the 
bellows. The bellows is a lightweight and flexible gastight 
seal between the driver and the tube. An aluminum plate is 
used as an interface between the reducer and the bellows to 
keep the bellows face rigid. For the STAR driver the reducer 
and plate are incorporated into a single unit machined from 
aluminum bar stock. 
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Figure II-3 shows a mechanical model of the 

electrodynamic transducer as a forced simple harmonic 
oscillator (SHO) , which consists of a moving mass (m.^) ; 
stiffness (k), mechanical damping or resistance (R m ) and 
driving force (F). Included in the moving mass are the 
masses of the voice coil, reducer, stiffener, and bellows. 
The stiffness includes the suspension stiffness of the voice 
coil surround, the stiffness of the gas volume trapped 
behind the diaphragm, and the stiffness of the bellows. The 
equations governing the SHO can thus be used to determine 
the moving mass and stiffness of the driver. The primary 
relationship used is : 

= k/m 0 , (11-19) 

or, the resonant angular frequency squared equals the 

stiffness divided by the moving mass. The use of this and 
other relationships will be discussed further in Chapter 

III, Driver Parameter Measurements. 

2 . Equivalent Circuit 

An electrical circuit can be drawn that is 
equivalent to an electromechanical system if the values 
given for R M , U, and C M in the circuit produce the same 
electrical effects as the electromechanical system itself. 
In the case of the electrodynamic driver, the equivalent 
circuit is shown in Figure II-4 with the values given by 
Kinsler, et al . (1982) as: 
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Figure II-3. Electrodynamic Transducer--Forced Simple 
Harmonic Oscillator Model. 
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Figure II-4. Equivalent Electrical Circuit for the 
Electrodynamic Transducer. 
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Rm = (Bl)=/R m , (II -20) 

L m = (Bl)=Vk, (11-21) 

C M = m 0 /(Bl)=, (11-22) 

where B1 is the transduction coefficient, R m is the 

mechanical resistance, and k and m, are the stiffness and 
moving mass as described previously. 

The properties of the equivalent circuit allow us to 
evaluate the performance of the electrodynamic driver as if 
it were an entirely electrical system. The amount of power 

the driver generates *_11 depend on the load it sees. Since 

one of the primary purposes of this thesis is to maximize 
the ratio of the acoutic power delivered to the electrical 
power provided (the electroacoustic efficiency), it is 

necessary to provide a load to the system in order to 
evaluate driver performance. This load is originally 
represented as a complex electrical impedence in the 
equivalent circuit. The actual value of the load will depend 
on the design of the resonator but can be characterized by 

an acoustic impedance which can be transformed to an 

equivalent electrical load impedance. 

3 . Losses 

In our search for the highest electroacoustic 
efficiency it is useful to understand where possible losses 
may occur and do what can be done to minimize them. Losses 
that are an intrinsic part of the electrical driver are eddy 
current losses in the magnetic structure generated by the 
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current in the voice coil. Joule heating losses due to the 
drive current dissipated by the voice coil electrical 
resistance, and mechanical losses due to the internal fric- 
tion of the surround and the viscosity of the gas in the 
"gap." These losses need to be taken into account. Very 
little can be done about the eddy current losses and the 
viscosity of the gas which is necessary to improve heat 
removal from the voice coil. The losses due to surround 
(suspension) mechanical resistance can be controlled by 
choice of the materials and removal of extraneous suspension 
material (see Section III.H). 

D. COUPLED SYSTEM 

1 . Computer Model 

Any discussion of the driver section of the 
thermoacoustic refrigerator is incomplete without some 
understanding of how the driver parameters are coupled to 
the resonator (acoustic load). The equivalent electro- 
mechanical circuit for the coupled system is the same as 
that shown in Figure II-4 with the acoustic load impedance: 

z i = -Pocsj [ (l + j<x/k)/(l+(«/k) 2 ) ♦ (11-23) 

( (cos(kL)sin(kL) + j sinh ( <xL ) cosh ( «L ) ) -r 

( sin= ( kL ) cosh 2 ( <xL ) + cos 2 ( kL ) sinh 3 ( «L ) ) ) ] , 
where Po i s the density of the gas, c is the speed of sound 

(labeled "a” in previous discussions), s is the cross- 

sectional area of the resonator, k is the wave number 
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(= u>/ c ) , « = u)o/(2Qc) , Q is the quality factor of the reso- 
nator, and L is the effective length of the resonator 
(Kinsler, et al . , 1982). This circuit combines the equiva- 
lent circuit of the electrodynamic transducer with the com- 
plex impedance load representing the resonator. One way to 
examine this coupling is through computer modeling. 

LT Michael P. Susalla, USN, developed a computer 
program which models the thermoacoustic refrigerator as the 

equivalent electrical circuit. The program models each of 
the mechanical and acoustical components as equivalent 
electrical components, as discussed above. Using the 
electrical equation from the combined circuit we determine 
the system performance as different parameters are varied. 

a. Input 

There are two sets of input parameters needed 
for the model: (1) driver parameters and (2) resonator 
parameters. The driver electrical parameters include the 
voice coil inductance (LE) and resistance (RE), and the 
current (I) supplied to the voice coil (which is held at one 
ampere for all test cases). The driver mechanical parameters 
used in the model are the same equivalent resistance, induc- 
tance, and capacitance as those developed for the equivalent 
circuit. Therefore the parameter values needed for input are 
the transduction coefficient (BL) , the mechanical resistance 
(RM), the moving mass (M) , the suspension stiffness (SSUP), 
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the volume of the gas behind the diaphragm (VOL), and the 
cross-sectional area of this volume (A) . 

The resonator parameters needed include the tube 
quality factor (Q) , the speed of sound (C), the resonance 
frequency (FO) , the density of the gas in the tube (RHO) , 
the resonator cross-sectional area <S) , and the equivalent 
length of the resonator (L) . For this model the piston area 
is set equal to the resonator cross-section (S) . For the 
tests reported here they were nearly equal but could differ 
sustantially in other designs, 
b. Output 

Four output plots are generated by the graph 
program that accompanies the model: (1) voltage across the 

driver vs. frequency (since we put unit current into the 
driver, this graph can also be read as the input electrical 
impedance), (2) piston velocity vs. frequency, (3) acoustic 
power delivered to the load (resonator) vs. frequency, and 
(4) electroacoustic efficiency vs. frequency. Each of these 
plots should produce peaks at both the driver and the tube 
resonant frequencies. 

The model uses the electrical equation for the 
coupled equivalent circuit to calculate voltage directly 
(holding the current constant at one ampere) for the first 
output plot. It also directly calculates the force on the 
moving mass, F = B1I. From this force, the piston velocity 
for the second plot is found by dividing the force by the 
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mechanical impedence (Z m »<=»-,). Pressure is calculated as tube 
impedence times volume velocity, and acoustic power is 
pressure times volume velocity times the cosine of their 
phase difference (Equation 1-3) for the third plot. For the 
final plot the electrical power is calculated as voltage 
times current times the cosine of their phase difference 
(Equation 1-7), and the model takes a ratio of the electri- 
cal to the acoustical powers for the electroacoustic effi- 
ciency. Since the components of these powers are complex and 
have phase differences, the actual powers are calculated as 
one half the real part of [current times complex conjugate 
of voltage] , and one half the real part of [pressure times 
the complex conjugate of volume velocity] . Each of the above 
mentioned plots is the discussed output parameter versus 
frequency . 

2 . Assumed Typical Parameters 

Table II-l lists the input parameters with sample 
values for the Altec 290-16K driver in air with a long 
straight resonant tube attached (see Chapter III). Sample 
outputs of the four plots using these input values are given 
in Figures II-5 through II-8. 

Figure II-5 plots the driver voltage as a function 
of frequency. Since the current is held constant at one 
ampere, this plot actually gives us information on driver 
electrical input impedance vs. frequency. Notice that there 
are two peaks in this curve. The first peak occurs at a 
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TABLE II-l. MODEL INPUT VALUES- -EXAMPLE DATA 



Symbol 


Meaning 


Units 


Value 


RM 


Mechanical resistance 


kg/s 


1 . 33 


M 


Moving mass 


kg 


0 . 0151 


SSUP 


Spring constant 


N/m 


9 . 2X103 


RHO 


Density of gas 


kg/m 3 


1.2 


C 


Sound speed 


m/s 


342 . 0 


s 


Resonator X-sec . area 


m 3 


9.35X10- 


BL 


Driver transduction factor 


N/A 


21 . 5 


Q 


Resonator quality factor 




60 . 0 


FO 


Resonator resonance freq. 


Hz 


240 . 0 


A 


Back volume X-sec . area 


m 2 


/ . 9 2 a 1 U 


VOL 


Back volume 


m 3 


1 . 0X10-= 


LE 


Voice coil inductance 


H 


1 . 0X10--* 


RE 


Voice coil DC resistance 


a 


10 . 6 


I 


Driver current 


A 


1 . 0 


L 


Resonator length 


m 


0 . 712 
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Figure II-5. Computer Model Sample Output--Driver 
Voltage vs. Frequency for One Ampere 
Current into Driver. 
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PISTON VELOCITY VS. FREQUENCY 




Figure II-6. Computer Model Sample Output--Piston 
Velocity vs. Frequency for One Ampere 
into Driver. 
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(xh) AQNunba 




Figure II-7. Computer Model Sample Outpu t--Acoustic 
Power vs. Frequency for One Ampere. 
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Figure I 1-8 . Computer Model Sample Output-- 
Electroacoustic Efficiency vs. 
Frequency . 
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lower frequency and has a much higher amplitude than the 
second. This higher amplitude peak occurs near the driver 
resonant frequency ( f 0 = ( 1/2* ) [ k/m„] 1 = 124 Hz). It 
dominates the curve because it is the driver impedance that 
we are plotting, and the voice coil motion, and consequently 
the "back emf , " is the largest at that frequency. The fact 
that there is a second peak near the tube resonance (FO = 
240 Hz) shows the effect of the coupling of the two parts of 
the system. 

In the plot of piston velocity vs. frequency (Figure 
I 1-6 ) , the dominant peak again occurs near the driver 
resonance for the reasons presented above. The smaller 
coupled peak at the tube resonance is due to the change in 
the load impedance which occurs near the tube resonance. 

Figure II-7 shows the plot of the acoustic power vs. 
frequency. Now the dominant peak occurs at the tube 
resonance because acoustic power is primarily a function of 
the tube parameters . Here the secondary peak at the driver 
resonance is a consequence of the greatly increased piston 
velocity at the driver resonance putting energy into the 
smaller non-resonant acoustic load of the tube. 

In the plot of electroacoustic efficiency vs. 
frequency, Figure II-8, the primary peak at the tube 
resonance is sharp and the secondary peak at the driver 
resonance is small but visible. It is important to note that 
the electroacoustic efficiency peaks at a value of nearly 
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60% even though the two individual resonances are separated 
by nearly an octave in frequency. The peak acoustic power 

for one amp as shown in Figure II-7 is 10.3 W. The losses 
(40%) are due primarily to Joule heating in the voice coil 
DC resistance ( P R = 1/2 (I 2 R DC ) = 5.3 W). This accounts for 
34%. ' The remaining loss can probably be attributed to 
mechanical resistance losses in the suspension. 

Appendix A contains a printout of the calculation 
program. Appendix B contains a printout of the graph 
program . 
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Ill . 



DRIVER PARAMETER MEASUREMENTS 



A. INTRODUCTION 

In Sections II. B and II. C we discussed the choice of the 
electrodynamic transducer to drive the thermoacoustic refri- 
gerator. Given this choice, we needed to develop a system 
for measuring the mechanical parameters of this type of 
driver to be used in determining the electroacoustic effi- 
ciency. The relevant parameters (as discussed in Section 
II. D) are: (1) moving mass, (2) suspension stiffness, (3) 
transduction coefficient (Bl), and (4) mechanical resis- 
tance. 

The driver we selected for the preliminary measurements 
is the Altec model 290-16K moving-coil loudspeaker. There 
are several pieces of equipment needed to interface the 
driver with the resonator. Three of these parts--the reducer 
cone, stiffener, and dynamic bellows--become an integral 
part of the driver, and their effects on the driver parame- 
ters must be taken into account. We first measured the 
parameters of the driver alone, then added the other parts 
and determined the new parameters for the entire system. The 
driver-bellows combination was then fitted to one end of a 
resonant tube. The actual resonator used in the thermo- 
acoustic refrigerator will be a complex structure that is 
not easily modeled. For the purposes of the driver parameter 
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measurements 



we used a resonant tube of uniform cross- 



section to allow easy comparison of experimental values to 
theory. The non-driver end of the closed tube housed several 
microphones. The electroacoustic efficiency of the system at 
resonance was determined from measurements of electric power 
delivered to the driver and the acoustic pressure at the 
microphone. The acoustic power dissipation was determined 
from the quality factor of the tube at resonance and the 
microphone sensitivity (the relationship between the pres- 
sure that the microphone senses and its voltage output). 



B. IMPEDANCE ANALYZER 

The primary instrument used in the testing and evalua- 
tion of the electrodynamic driver was the Hewlett Packard 
Model 4192A LF Impedance Analyzer. According to the manufac- 
turers Operation and Service Manual, the 4192A is a fully 
automatic, high performance test instrument designed to 
measure a wide range of impedance parameters as well as 
gain. phase, and group delay. Its two measurement display 
sections provide direct readout of the selected measurement 
parameters with 4.5 digit resolution along with the appro- 
priate units. In normal mode operation, the 4192A performs 
approximately five measurements per second. It also provides 
an average measurement mode (one measurement per second) for 
higher resolution and noise rejection. A high speed mode can 
perform ten measurements per second. 
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The 4192A can provide measuring frequency, OSC voltage 
level, and DC bias voltage (impedance measurements only) 
equivalent to actual operating conditions . The sweep capa- 
bility of the built-in frequency synthesizer and DC bias 
source permits quick, convenient, and accurate measurements. 
Measuring frequency or DC bias voltage can be automatically 
or manually swept in either direction. OSC level can be 
manually swept in either direction in lmV increments. Actual 
test voltage across or test signal current through the 
device under test is also measured . 

In amplitude/phase measurements , the 4192A can measure 
four transmission parameters - gain/loss. level, phase, and 
group delay. In impedance measurements, the 4192A can 
measure eleven impedance parameters - absolute value of 
impedance, absolute value of admittance, phase angle, resis- 
tance, reactance, conductance, susceptance, inductance, 
capacitance, dissipation factor, and quality factor. All 
impedence measurements have a basic accuracy of 0.1% and a 
resolution of 4.5 digits. Also, the circuitry of the 4192A 
provides direct and accurate impedance measurements of both 
grounded and floated devices. 

C. MECHANICAL PARAMETERS 

The most important parameters for a moving coil trans- 
ducer are the moving mass, suspension stiffness, and the B1 
(transduction) coefficient. The moving mass of the driver 
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includes the mass of the voice coil, the diaphragm, and the 
hydrodynamic mass of the air being pushed by the diaphragm. 
The masses of the reducer and stiffener are also important, 
as are the moving mass and suspension stiffness of the 
bellows. The moving mass and stiffness of the combined 
driver-reducer-bellows system is also important. The values 
for the combined system should equal the sum of the values 
of the individual parts. The stiffness, mass, and Bl coeffi- 
cient give us the information we need to construct the 
equivalent electrical circuit, as discussed in Subsection 
II. C. 2. Also, we have to insure that the stiffness of the 
driver system is just large enough to produce a mechanical 
resonance at the tube resonant frequency, at which point the 
driver piston looks like a soft (or open) tube end. This may 
appear to be contradictory to having a velocity node of the 
standing wave at the piston, but it is not. It is how we get 
high efficiency. 

The important design parameters of the tube are its 
resonance frequency and quality factor. The length of the 
tube and the speed of sound in the contained gas determine 
the resonance frequency. 

1 . Dynamic Mass Loading of Driver 

We used the driver’s similarity to a simple harmonic 
oscillator to find the moving mass and suspension stiffness 
of the driver by taking advantage of the fact that the 
resonance frequency squared equals the stiffness divided by 
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the mass . 



If we add mass to the driver, it should linearly 



change the inverse of the frequency (or period) squared: 

T= = l/f= = 4Tt=(m 0 + m)/k, (III-i) 

where m^ is the moving mass, m is the added mass, and k is 

the suspension stiffness. The instrument we used for this 
measurement was the HP 41Q2A Low Frequency Impedance 
Analyzer, described in section B. We added various masses of 
putty (one to nine grams) to the diaphragm and found the 
resonance frequency of the driver for each mass. The putty 
mass was measured on the Sartorius model 2403 mass balance 
(serial number 1287). This balance has an accuracy of ± 
0.00005 g. However, the putty could be expected to gain or 
lose more mass than that due to handling. Based on mass 
measurements of the putty after use, we expect an 
experimental accuracy of ± 0.005 g for the masses. A more 
significant uncertainty involving the masses arises from the 
inability to place the masses evenly on the voice coil or to 
produce a perfectly rigid coupling between the masses and 
the voice coil. It is not possible to assign a specific 
numerical percentage error to this uncertainty, but we must 
take it into account when evaluating the data associated 
with mass loading. A typical plot of the period squared vs. 
added mass in air is shown in Figure III-l. We used the 
least squares method to fit a straight line through the 
data. We determined the suspension stiffness (k = 4;rt=/ s ]_ 0 p e ) 
to be : 
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PERIOD SQUARED (x10-6 S*S) 



ALTEC 290 - 1 6 K MASS LOADING 




ADDED MASS (g) 



Figure III-l. Least Squares Fit of Period Squared vs. 

Added Mass to Determine Stiffness and 
Moving Mass of the Altec 290-16K Driver. 
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k 



2.7 



0 . 1 X 10-* N/m. 



( 111 - 2 ) 



and the moving mass (m Q = intercept/slope) : 

m 0 = 5.8 ± 0. 1 g, ( I I I -3 ) 
for this driver in air. 

The moving mass of the driver in vacuum is expected 
to be less than the value in air because in vacuum the 
driver is not pushing a hydrodynamic mass. However, we did 
not take any measurements of this driver in vacuum. 

2 . Bellows Parameters 

The Servometer custom-designed dynamic bellows 
(shown in Figure III-2 with test stand) is used to provide a 
flexible gas seal between the driver and the resonator. 

The stiffness of the bellows was measured using two 
independent methods - static and dynamic. For the static 
method we used Hooke’s law: k = -F/x. The force was deter- 
mined by pulling a partial vacuum on the bellows and multi- 
plying the vacuum pressure (read on a manometer) by the 
surface area of the bellows (A = 1.1 ± o.i X 10~ 3 m = ) We 
used a G. L. Collins Corporation model LMS-1 19V22-01 Linear 
Variable Differential Transformer (LVDT) (serial 274674) 
connected to an EG&G 128A Lock-in Analyzer to determine the 
displacement of the bellows under a partial vacuum. Figure 
III-3 shows the LVDT in its calibration stand. The calibra- 
tion curve for the LVDT is shown in Figure III-4. The least 
squares fit slope of this curve gives us the calibration 
constant for the LVDT : 
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Figure III-2. Photograph of Servoraeter Custom-Designed 

Dynamic Bellows — With and Without 
Test Stand. 
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Figure III-3. Photograph of Linear Variable 

Differential Transformer ( LVDT ) 
in its Calibration Stand. 
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LOCK-IN OUTPUT VOLTAGE (mV dc) 



LVDT CALIBRATION 




Figure III-4. Calibration Data for LVDT Voltage 

Output vs. Displacement. 
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BELLOWS DYNAMIC MASS LOADING 




Figure III-5. Least Squares Fit of Period Squared 

vs. Added Mass for Dynamic Loading 
of Bellows. 
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dV/dx = -0.378 ± 0.0005 V/in = -14.9 ± 0.02 V/m. (III-4) 
Table III-l gives the static test data for the stiffness 
measurement. By averaging the values we obtained: 

k = 4.3 ± 0.2 X 10° N/m. (II 1-5) 

TABLE III-l. BELLOWS STIFFNESS — STATIC 
MEASUREMENT RESULTS 



DT Voltage 


Pressure 


Pressure 


Force (=P*a) 


Stif f ness 


(mV) 


(cm Hg) 


(XlO^Pa) 


(N) 


(X10- 3 N/m) 


14 . 84 


2 . 9 


3.87 


4.39 


4.40 


15.62 


3.1 


4 . 13 


4 . 69 


4 . 47 


12.89 


2 . 3 


3.07 


3.48 


4 . 02 



The dynamic measurement was a mass loading technique 
similar to that used for the driver. In this case, we 
measured the period as a function of added mass using an 
Endevco model 2225M2 accelerometer (serial FA80) and a 
Nicolet model 3091 digital storage oscilloscope to record 
the free decay of the bellows after it had been disturbed by 
the blunt wooden end of a cotton swab. We then plotted 
period squared versus added mass and found the stiffness and 
moving mass from the slope and intercept of the least 
squares fitted line (see Figure III-5) as discussed in the 
previous Section. This gave us values of: 
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k 



(II 1-6 ) 



= 4.56 ± 0.05 X 10 3 N/ra 
and » 0 = 5.84 ± o . 05 g. (III-7) 

This value for moving mass includes the accelerometer cable 
and is not representative of the actual moving mass of the 
bellows. Subtracting the weight of the accelerometer gives 
us a value of approximately: 

m 0 = 1.5 ± 0.1 g. (Ill -8 ) 

for the bellows. We know the moving mass has to be less than 
2.3 g, since that is the bellows' total mass— moving and 
static. The values obtained for bellows stiffness using the 
two different methods were within experimental error. We 
therefore felt confident using a value of k B = 4.5 ± 0.1 X 
10 3 N/m in our later calculations. 

3 . Combined System 

After we glued the reducer (m = 4.8 g) . stiffener 
plate (m = 4.4 g) and bellows to the driver, we wanted to 
find the stiffness and moving mass of the combined driver- 
reducer-bellows system. We first tried a dynamic mass 
loading measurement (same procedure as with the driver 
alone), but found that the period squared versus added mass 
plot (see Figure III-6) gave unrealistic results. So instead 
we used the relationship: 

F = -kx = B1I (II I -9 ) 

or k = -BlI/x = B1 (di/dV) (dV/dx) (III-10) 

and made a static measurement using the LVDT. We expected 
the additional equipment (reducer, stiffener, and bellows) 
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MASS LOADING DRIVER-CONE-BELLOWS 




ADDED MASS (g) 

Figure III-6. Least Squares Fit of Period Squared 

vs. Added Mass for Combined Driver- 

Reducer-Bellows System. 
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ALTEC 290-16K STATIC DISPLACEMENT 




CURRENT IN VOICE COIL(mAdc) 

Figure III-7. Driver Displacement vs. Voice Coil 

Current as Measured with LVDT/ 
Lock-ir. Analyzer. 
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to change the moving mass and suspension stiffness of the 
system but not the transduction coefficient ( B1 ) . We applied 
a DC current to the driver and recorded the LVDT voltage 
output . then changed the current and recorded the new vol- 
tage. This process produced the values: £i = 283 ± 0.5 mA 

and &v = 2.58 ± 0.005 mV. We substitute these values and the 
value of B1 = 21.5 ± 0.5 N/A (experimental measurement of 
this quantity is discussed in the next section) into Equa- 
tion III-10 for a combined stiffness value of: 

k = 3.5 ± 0.1 X 10 4 N/m. (III-ll) 

If we add the stiffness values for the driver and bellows 
together a total stiffness of k = 3.2 + q.i y N / m r,-, r 

the system is obtained. Thus the total experimental stiff- 
ness for the combined system differs by about 10% from the 
sum of the individual stiffnesses. Using the relationship 
= k/(2rtf) = . and a measured resonance frequency of 

235 ± 0.5 Hz, the combined moving mass is: 

me = 16.2 ± 0.5 g, 

compared to m 0 = 15.9 ± 0.5 g for the 
added together . The two types of mass 
within 2%. 

D. TRANSDUCTION COEFFICIENT (BL) 

We determined the B1 coefficient of the driver using 
three independent methods. The first method was to remove 
the voice coil from the magnet and directly measure the 



( III-12) 

individual masses 
measurement agree 
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magnetic induction using a Dowdy RFL Industries model 912 
Guassmeter (serial 808) using a Hall effect probe model 
912015 (serial 21122) (B = 2.2 ± 0.05 Tesla). We then 
measured 1 (the length of the wire) by counting the number 
of turns of wire in the coil (under a microscope) and 
multiplying this number (44) by the circumference of the 
coil (1 = 44 X rt X diameter = 9.98 ± 0.005 m) . Combining 
these two measurements gave us the value: 

B1 = 22.0 ± 0.7 N/A. (III-13) 
In the second method we used the relationship: 

F = B1I = -kx or B1 = -kx/I = -k (dx/dV) (dV/dl ) . (III-14) 
We previously determined k for the driver alone (Equation 
III-2) , and now used the LVDT and lock-in analyzer to find 
dx/di. The LVDT calibration was discussed in the previous 
Section, with a calculated value of dV/dx = -14.9±0.05 V/m. 
We then measured the diaphragm displacement as a function of 
DC-current through the voice coil. This gave us a value of 
dV/dl = 11.75 + o.05 mV/A (see Figure III-7). Using the 
relation in Equation III-14 we obtain the value: 

B1 = 21.2 ± 0.1 N/A. ( I I I -15 ) 
In the third method we used the relationship: 

F = Bli = mg or B1 = mg/i. (III-16) 
For this we measured the change in DC-current through the 
coil as a function of added mass, using the LVDT as a 
monitor to keep the displacement of the diaphragm constant. 
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Three different mass/current combindtions ( Tabl*=: III-2) 
produced an average value of : 

B1 = 21.4 ± 0.5 N/A. (I II -17) 

Averaging the results of these three methods gave us a 
value of B1 = 21.5 ± 0.5 N/A with 2% deviation among the 
three methods . 

TABLE III-2. STATIC MEASUREMENT OF ALTEC 290-16K 
TRANSDUCTION COEFFICIENT (BL) 



Am ( g) 


AI (mA) 


B1 = mg/I 


132 . 5 


64 . 4 


20 . 2 


157.5 


74 . 0 


20 . 9 


257.5 


108.2 


23 . 3 



E. MICROPHONE CALIBRATION 

Microphone calibration is an important part of the 
process leading to the measurement of the electroacoustic 
efficiency. A microphone is used to measure the acoustic 
pressure inside the resonator by converting the acoustic 
pressure to an output voltage. The ratio of the open circuit 
voltage to the pressure is called the open circuit 

microphone sensitivity. 

We used three microphones for the calibration.- (1) a 
General Radio (GR) model 1962-9611 half-inch electret 
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condenser microphone (serial 11811) with GR preamplifier 
model 1560-P42 (serial 4061); (2) an Endevco model 8510B-5 
piezoresistive bridge microphone (serial 73MJ); and (3) an 
1/8 inch B&K model 4138 microphone (serial 19683) . We first 
calibrated the GR using a GR type 1562-A sound-level cali- 
brator (serial 19683, 114dB SPL , 250 Hz), and obtained a GR 
sensitivity of = 7.42 mV/Pa. This compares within about 
1 % of the calibration certificate provided by the manufac- 
turer of 7.50 mV/Pa. Additional checks were made with a B&K 
pistonphone calibrator which substantiated the GR calibra- 
tion results . 

To calibrate the other two microphones we attached the 
resonator tube to the 290-16K driver. We used an HP 3325A 
Synthesizer/Function Generator to excite the driver at 240 
Hz, the resonant frequency of the tube. We recorded the 
output voltages for all three microphones and used the 
comparison relationship: 

M i/V x = M = /v = (III - 18 ) 

to determine the sensitivities of the B&K and Endevco. The 
B&K sensitivity was determined to be Mo = 0.463 mV/Pa, 
comparing within 2% of the manufacturer's specification 
value of 0.470 mV/Pa (adjusted to include the capacitance of 
the B&K model VA0160 1/8” to 1/4" adapter with an expected 4 
dB "line loss"). The sensitivity of the Endevco microphone 
was calculated to be 8.17 X 10-* v/ Pa, comparing within 1% 
of the manufacturer's specification of 8.09 X 10-.& y /Pa For 
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all three microphones we measured sensitivity values within 
2% of the manufacturers’ specifications. 

F. ELECTROACOUSTIC EFFICIENCY 

The electroacoustic efficiency is the ratio of the 
acoustic power dissipated in the resonator to the electric 
power delivered to the driver. The electric power was 
determined using the relationship: 

P-x«= = V I cos S, ( II I -19 ) 
where V is the voltage across the driver terminals, I is the 
current flowing in the driver, and $ is the phase angle 
between the voltage and the current. The values for voltage, 
current and phase angle delivered to the driver were 
obtained using the HP 4192A. 

Calculation of the acoustic power was a little more 
complex. We started with the basic relationship for the 
quality factor of the tube: 

Q = 27tE^x/Ei_c , ( 1 1 1-20 ) 
where E ST i s the energy stored in the tube and Ei_ c is the 
energy lost per cycle. The acoustic power generated in the 
tube is the energy lost per cycle times the number of cycles 
per second, or frequency: 

P^c = Ei_ c f = 2nf Es-r/Q . (1 11-21) 
Following Burmaster (1985) the energy stored in a plane 

standing sound wave is equal to the energy density 
integrated over the volume of gas in the tube (VOL) : 
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(P 



/2;iP 0 ) * VOL . 



( III-22) 



Est = 

where P^ m ^ is the root mean square acoustic pressure in the 
tube, v is th e ratio of the specific heats of the gas in the 
tube (1.4 for air), and P 0 is the ambient pressure (one 
atmosphere or 1.013 X 10 3 Pa for these measurements). The 
r.m.s. acoustic pressure is the r.m.s. microphone output 
voltage divided its sensitivity (V/M) as discussed in the 
previous Section. The expression for the acoustic power 
becomes : 

P«c - <rtfV=VOL)/(QM=*Po) . (Ill -23 ) 

The only remaining quantity in Equation III-23 that 

needs to be calculated separately from the efficiency 
measurement is the quality factor of the tube (Q) . To 
determine Q we removed the 290-16K driver from the resonator 
and replaced it with an electret transducer to drive the 
tube. This was necessary because there were losses due to 
the 290-16K that made the tube Q appear to be lower than it 
actually is (we discovered this problem while investgating 
the reasons for the calculated efficiency to be greater than 
one). We used the B&K microphone for the pressure (converted 
to voltage) readings at the other end of the tube. We plot- 
ted voltage as a function of frequency, and determined the 
resonant frequency (f Q ) and frequencies of the half power 
points (f^ and fh) . Using the relationship: 

0 = f«,/(f„-f x ) ( III-24 ) 

we obtained a value of Q = 60 ± 1 for the tube. 
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To determine the electroacoustic efficiency we returned 
the tube to the 290-16K driver (see Figure III-8) and drove 
the tube at resonance. We recorded the voltage, current, and 
phase angle for the driver and the B&K output voltage. The 
electroacoustic efficiency was then calculated using: 

n E « = Pac/P-x-c. (II 1-25 ) 
Table III-3 lists the electric power, acoustic power, and 
efficiency for five test runs at driver currents ranging 
from 3-38 mA (with corresponding voltages 0.15-1.43 V). The 
electroacoustic efficiency ranged from 84.2 for the lowest 
driver current to 75.9 for the highest. 

TABLE I I 1-3 . ALTEC 290-16K ELECTROACOUSTIC 
EFFICIENCY MEASUREMENTS 



P-i~= (mW) 


P<c*o ( mW ) 


Hea (%) 


0.56 


0 . 47 


84.2 


3.41 


2.70 


79.2 


10 . 0 


7.79 


77.9 


17 . 0 


13. 1 


77.2 


54.2 


41 . 2 


75.9 



G. COMPARISON 


TO COMPUTER MODEL 




The final 


task in the 


evaluation 


procedure for the 


electrodynamic 


driver is to 


compare the 


measured electro- 



acoustic efficiency to that predicted by the computer model 
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Figure III-8. Photograph of Altec 290-16K Driver 

Assembly with Resonator Tube and 
Endevco Microphone. 
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discussed in Section II. D. We have already measured or 
calculated all the input values for the model except for the 
driver’s mechanical resistance (R m ). 



1 . Mechanical Resistance 

We can determine R m by using the relationship: 

Rm = 2rtf«dWQ, ( 1 1 1 -26 ) 

where f Q is the resonant frequency, nu= is the moving mass 
and Q is the quality factor of the driver. There are several 
methods available for determining the driver’s quality 
factor. One is to plot the magnitude of the modified driver 
impedance vs. frequency and pick the half power points off 
the graph. Then, as with the tube, 0 is given by: 



Q = f 0 /(f„-fx ) , (III -27 ) 
where f and f i are the frequencies associated with the 
upper and lower half power points (located at [ Z yn, 
where [Z] is the magnitude of the driver impedance which is 
modified by subtracting the direct current (DC) resistance 
from the measured complex impadance) . 

Another method is to use the relationship [Hofler, 



1986] : 



Q = £• fo dO/df, (III- 28 ) 
where is the phase angle between driver resistance (R) and 
reactance (X), the two components of modified impedance (Z = 
R + iX . where R is the measured resistance minus the DC 
(direct current) resistance). This modified impedance is 
proportional to the voice coil velocity. The quantity d$/df 
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is the slope of the phase angle versus frequency curve 
evaluated at the resonant frequency (where the curve is a 
straight line) . We used the HP4192A to obtain the data for 
these measurements . 

Figure III-9 shows a plot of modified impedance vs. 
frequency for the Altec 290-16K in air. From this graph we 

find f o=235 . 2 Hz, f^=249.6 Hz and f» =222.0 Hz. The Q for 
these values is 8.52. Figure III-10 shows a plot of phase 
angle vs. frequency for the same data set. Over this large 
range of frequencies the phase angle curve approximates an 
arctangent. At resonance, however, where ds/df is evaluated, 
the curve is linear (see Figure III-ll). Equation III-28 is 
the slope of the arctangent at resonance (f Q ) . The slope of 
the plotted data is 71.8X10 -3 sec (dO was converted from 
degrees to radians). The Q calculated from this method is 
8.44. These values for Q agree within 1%. From these values 
we get R m =2 . 77 kg/s. 

We tried yet another method for determining driver 
Q,. This third method measures the decay constant, x, using a 
digital storage oscilloscope to record the "ring down" of 
the driver. This technique takes less time to run the exper- 
iment and reduce the data using the relationship: 

Q = mX/2 or X = 2Q/a> = QT/tf (III-29) 
where <j> is the resonant angular frequency in radians/sec, T 
is the period in seconds, and t is the free decay time 
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ALTEC 290-1 6K IMPEDANCE 




Figure III-9. Altec 290-16K Modified Impedance 

Magnitude vs. Frequency. 
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Figure III-10. Altec 290-16K Modified Impedance 

Phase Angle vs. Frequency--Full 
Frequency Range. 
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Figure III-ll. Least Squares Fit of Altec 290-16K 

Modified Impedance Phase Angle as a 
Function of Frequency Near Resonance. 
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constant in seconds. This method will be discussed in 
greater detail in Chapter IV. 

2 . Model Results 

Table III-4 lists the computer model input values 
for the Altec 290-16K. Figures III-12 through III-15 show 
the output plots for these values as discussed in Section 
II. D. In this case the resonance peaks are much closer 
together, and the magnitude of the secondary peak is closer 
in value to the primary peak. Once again the driver reso- 
nance is at a lower frequency than the resonator resonance. 
It is interesting to note that the modeled efficiency (about 
60%) is 15-20% less than the measured efficiency reported in 
the previous Section. We believe that this results from the 
fact that for the computer model, a high tube O means that 
the driver losses, represented by R m , play a more important 
role in the efficiency calculation, and the results are very 
sensitive to small errors in R m . The primary ways for the 
driver energy to be dissipated are: (1) Joule heating due to 

DC resistance in voice coil and eddy currents induced in the 
magnet structure due to the time-varying magnetic fields 
generated by the AC voice coil current (small but not negli- 
gible) ; (2) mechanical resistance of driver (R m ); and (3) 

energy transferred to acoustic load. Another possible reason 
for the different efficiencies is that the efficiency is not 
independent of current hence the real physical system has 
some non-linearity, as Table III-3 clearly shows. The 



79 



TABLE III-4. COMPUTER MODEL INPUT --ALTEC 290-16K DRIVER 



Symbol 


Meaning 


Units 


Altec 


RM 


Mechanical resistance 


kg/s 


2.77 


M 


Moving mass 


kg 


0 .0162 


SSUP 


Spring constant 


N/m 


3 . 52X1 0 - * 


RHO 


Density of gas 


kg/m 3 


1 . 2 


C 


Sound speed 


m/s 


342 . 0 


S 


Resonator X-sec . area 


m = 


9 . 35X10"“* 


BL 


Driver trans. factor 


N/A 


21 .5 


0 


Resonator quality factor 




60 . 0 


FO 


Resonator resonance freq. 


Hz 


240.0 


A 


Back volume X-sec . area 


m = 


4 . 13X1 0~ 3 


VOL 


Back volume 


m 3 


2 . 3X10-* 


LE 


Voice coil inductance 


H 


1 . oxio- 3 


RE 


Voice coil DC resistance 


ft 


10 . 6 


I 


Driver current 


A 


1 . 0 


L 


Resonator length 


m 


0 . 712 
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Figure III-12. Computer Model Output for 

Altec 290-16K--Driver Voltage 
vs. Frequency for One Amp. 



81 





Figure I I I - 13 . Computer Model Output for 

Altec 290-16K--Piston Velocity 
vs. Frequency for One Amp. 



82 



ACOUSTIC POWER VS. FREQUENCY 




Figure III-14. Computer Model Output for 

Altec 290-16K — Acoustic Power 
vs. Frequency for One Amp. 
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ELECTROACOUSTIC EFFICIENCY VS. FREQUENCY 




Figure III-15. Computer Model Output for 

Altec 290- 16K--Electroacoustic 
Efficiency vs. Frequency. 
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computer model is based on a linear theory and will have an 
efficiency which is independent of drive current. 

H. DRIVER RESONANCE ALTERATION 

The refrigerator must operate at the resonant frequency 
of the resonator to obtain the highest electroacoustic effi- 
ciency and to maintain stack position with respect to the 
standing acoustic wave (discussed in Sub-section II. A. 4). It 
appears from the computer model results and our experiments 
so far that the operating frequency should also be near the 
driver's resonance to get the most power out of the driver. 
The resonant frequency of the resonator is determined by the 
effective length of the resonator and the speed of sound in 
the contained gas. The effective length for a tube of uni- 
form cross-section and closed end is the actual measured 
length. For the STAR resonator, however, the effective 
length is the length that would produce the measured reso- 
nant frequency for the given sound speed if the resonator 
were a uniform cross-section closed-ended tube (see Figure 
1-2 for approximate shape of resonator). The sound speed can 
be altered by changing the percentage of xenon in the 
helium-xenon mixture. since the sound speed depends on the 
density of the mixture, and the density depends on the 
percentages of the individual gases present (see Equations 
IV-24 and IV-26). The gas mixture and effective resonator 
length will be chosen to maximize the coefficient of 
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performance of the refrigerator. Therefore if we want to 
match the resonances, we will have to alter the resonance of 
the driver once the resonant frequency of the resonator has 
been determined. 

The resonant frequency of the driver is proportional to 
the square root of the stiffness divided by the moving mass. 
If we want to raise the resonant frequency of the driver we 
can either increase the stiffness or decrease the mass. 
Conversely, to lower the resonant frequency we can either 
decrease the stiffness or increase the mass. Generally we 
prefer to decrease stiffness and mass to get a larger 
velocity. Reduction of stiffness usually also reduces the 
mechanical resistance. 

For the Altec 290-16K attached to the uniform cross- 
section resonant tube in these preliminary experiments, the 
frequencies were already well matched (note the closeness of 
the peaks in Figures III-12 through III-14) and changing the 
driver frequency was not required. However, we know that in 
future applications we may need to alter, most probably 
lower, the driver resonant frequency. We developed an 
experimental procedure to do this . 

Adding moving mass to the driver would appear to be the 
easiest way to lower the driver resonance. Unfortunately the 
added mass would also reduce the velocity by increasing the 
mecahnical impedance above resonance. We therefore chose to 
experiment with decreasing the driver stiffness. Recall that 
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the driver-reducer-bellows system stiffness includes the 
stif f f ness of the bellows (4.5 X 10 13 N/m) and the suspension 
stiffness of the driver phenolic surround (29 X 10-" 3 N/m). 
The driver suspension was the obvious choice for 
experimental alteration for two reasons: (1) it is over six 
times as large as the bellows stiffness, and (2) it would be 
very difficult for us to change the bellows stiffness 
without affecting its performance. 

Since the 290-16K voice coil surround is made of a 
phenolic material, it was relatively easy to cut through it 
with an ordinary scalpel. We made a series of successive 
cuts in the surround material and ran measurements of the 
driver resonant frequency after each cut. The results of 
this experiment are given in Table III-5. We plotted the 
phase angle between the driver's electromotive force (emf) 
and the voice coil displacement as a function of frequency. 
For the emf we measured the voltage across the driver termi- 
nals in channel A of the HP 4192A Impedance Analyzer and 
assumed that the emf was in phase with the driver current 
and hence force. We used the photonic sensor to measure the 
displacement by connecting its voltage output to channel B. 
The 4192A then displayed the phase angle between the two 
voltages representing force and displacement. The driver is 
at resonance when this phase angle equals 90 degrees . The 
quality factor (Q) is equal to one-half the resonant fre- 
quency times the slope of the phase angle vs. frequency 
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curve evaluated at resonance (see Equation III-28). Using 
the previously determined values for moving mass (16 g) and 
bellows stiffness, we calculated the suspension stiffness: 

^ w — h fcofc kb«l lows ( I 11-30 ) 

where . We calculated the mechanical 
resistance (R m ) using equation III-26. We stopped cutting 

the surround when its stiffness was approximately equal to 
that of the bellows and therefore the total stiffness was no 
longer suspension dominated. 



TABLE III-S. RESULTS OF 
DRIVER RESONANT FREQUENCY ALTERATION 



Percent Surround 


Resonant 


Quality 


Suspension 


Mechanical 


Remaining 


Frequency 


Factor 


Stiffness 


Resistance 


( approx . ) 


(Hz) 




( XlO^N/m ) 


(kg/s ) 


100 


241 . 0 


8.55 


32 . 2 


2.33 


100 (six slits) 


229.5 


8 . 09 


28 . 8 


2 . 85 


80 


213.5 


5 . 17 


24 . 3 


4 . 15 


50 


183.5 


8 . 65 


16 . 8 


2. 13 


20 


119.5 


8.56 


4.52 


1 . 40 


Despite the 


anomoly at 


the 80% 


data point , 


we feel tha' 


this experiment 


shows that 


partial 


removal of 


the surroum 



can reduce stiffness, resonant frequency, and mechanical 
resistance. It appears from the above data that the flexing 
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of the phenolic surround is responsible for the majority of 
the mechanical losses in the driver. 
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IV. 



star driver parameters 



A. INTRODUCTION 

This chapter discusses the selection of the electro- 
dynamic driver, describes the parameter measurements on the 
chosen driver, and estimates the electroacoustic efficiency 
of the STAR using the computer model described in Section 
II . D. 

1 . Interface with NASA GAS Cannister 

In addition to the maximum weight of 200 pounds 
discussed in Section I.B, restrictions due to the Get Away 
Special (GAS) cannister dimensions give us a maximum payload 
height (28.25 inches) and diameter (19.75 inches). In our 
alloweu height we must fit two battery layers, the driver 
with its housing, the resonator assembly, and its associated 
vacuum space and thermal insulation. The battery layers are 
each about four inches high, for a total battery height of 
eight inches . This leaves about 20 inches for the driver and 
resonator assemblies (see Figure 1-4 for approximate 
configuration). The electronics package and gas reservoir 
should fit in the area radially surrounding the resonator, 
and will therefore not add to the height. Although the 
diameter of the cannister volume is 19.75 inches, the driver 
housing must attach to the standard GAS can lid at a 9.5 
inch diameter bolt circle, otherwise we would need a lid 
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specially designed for our project. As far as weight is 
concerned, approximately every pound we save elsewhere 
allows us to carry an extra battery cell, which results in 
an extra ten watt-hours of available energy. 

2 . Maximizing Electroacoustic Efficiency 

In addition to the space and weight restrictions, we 
want to maximize the electroacoustic efficiency of the 
driver. We can do this by obtaining the maximum value for 
the combination of parameters Bl/R m m 0 . A larger B1 will 
give us a larger driving force (F=Bli), which in turn gives 
us a larger acoustic volume velocity, one of the components 
of acoustic power. For a given value of magnetic induction 
(B) , a longer voice coil (achieved by having a larger voice 
coil diameter with the same number of turns. or the same 
diameter and more turns) will result in a larger B1 . This 
longer 1 would also result in a greater voice coil 
electrical resistance, so the two effects cancel. We also 
want to minimize the moving mass (m,.,) to allow greater 
velocities and minimize the driver's mechanical resistance 
(R«n) in order to minimize mechanical losses in the driver. 

3 . Driver Selection 

We had two off-the-shelf drivers available for con- 
sideration. The Altec 260-16K was discussed in Chapter III. 
The other driver was a JBL 2445J. Both drivers were designed 
for applications somewhat different from the STAR. Figure 
IV-l shows a side view of the 2445J driver. The hatched 
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Figure IV- 1 . Schematic Drawing of JBL 2445J Driver. 

Hatched Area Shows Unnecessary Throat 
Section . 
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area, known as the throat of the driver, is considered 
unnecessary for our application, and adds a lot of extra 
weight to the driver. The Altec 290-16K has a similar excess 
throat volume. We contacted JBL and asked if they would 
custom-design a driver which would eliminate this excess 
weight. Fortunately for our project, JBL was in the process 
of designing a new driver, the 2450 J, which uses neodymium- 
iron-boron magnets instead of the traditional ceramic. They 
agreed to custom-design a version of the 2450 J for our use. 
The resulting driver (hereafter referred to as NIB) is shown 
in Figure IV-2. The NIB driver accepts the same voice coils 
(four inch diameter, aluminum wire, and titanium suspen- 
sion/diaphragm) as the 2445J (ceramic magnet) driver. JBL 
claims that it should also have the same properties , 
including B1 . 

In the remainder of this chapter we describe the 
parameter tests on the NIB driver. We then compare the NIB 
driver to the 290-16K driver. The results of the 2445J 
parameter measurements are given for comparison. Finally we 
present results generated by the computer model for the NIB 
driver given various acoustic load impedances. 

B. NIB DRIVER PARAMETERS 

For the measurement of the NIB driver parameters, we 
used the techniques we developed in Chapter III for the 
measurement of the 290-16K driver parameters. 
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Figure IV-2. Photographs of JBL Neodymium- 
Iron-Boron Driver. 
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1 . Moving Mass and Suspension Stiffness 

The moving mass and suspension stiffness were 
calculated using the mass loading technique. We placed small 
masses (0.5 to 3.5 grams) of black putty on top of the 
voice coil in half -gram increments. For each added mass the 
resonant frequency of the driver was found using the HP 
4192A LP Impedance Analyzer (see Section III.B for a 
description of this instrument) . Figure IV-3 shows a sample 
plot of period squared vs. added mass for the NIB driver in 
air. The stiffness was found by calculating four pi squared 
divided by the slope. The moving mass is equal to the 
intercept divided by the slope. For the NIB in air these 
values are: 

k = 9.8 ± 0.1 X 10-* kg/s 2 ( IV-1 ) 

and m 0 = 7.8 ± 0.1 g, (IV-2) 

with a resonant frequency of 564 Hz. We repeated the 
measurement in a vacuum chamber and found: 

k = 8.6 ± o.l X 10^ kg/s 2 (IV-3) 

and m 0 = 3.2 ± 0.1 g, (IV-4) 

with a resonant frequency of 833 Hz (see Figure IV-4). 

The moving mass is greater in air because of the 
driver having to move the air itself. The discrepancy in the 
stiffness is not as easy to explain. To determine the actual 
value of k, we used the static mass loading method. This 
gave us a value of : 
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PERIOD SQUARED (S*S) 



NIB MASS LOADING • AIR 




Figure IV-3. Least Squares Fit of Period Squared vs. 

Added Mass for NIB Driver in Air. 
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Figure IV-4 . Least Squares Fit of Period Squared vs. 

Added Mass for NIB Driver in Vacuum. 
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6.4 ± 0.2 X 10"* kg/s = . 



k 



( I V-5 ) 



(The static loading measurements will be discussed in the 
next Sub-section.) 

This last value for suspension stiffness differs by 
about 40% from the first two. We did not have sufficient 
time to test this driver until we were sure we had the 
"correct" value for the stiffness. We recommend that the 
follow-on students do so. Lacking the correct value, we ran 
the computer model for two sets of data: one set represen- 

ting the "high" k (9.0 X 10"* kg/s 2 ) and one using the "low” 
k (6.4 X 10-* kg/s 2 ) . 

2. Transduction Coefficient ( B1 ) 

The transduction coefficient (Bl) was calculated 
using two methods. The first method was to measure B (magne- 
tic induction in the voice coil gap) and 1 (length of the 
voice coil) directly. We measured B with the Dowty RFL 
Industries model 912 Guassmeter (serial 303) using a Hall 
effect probe model 912015 (serial 21122) . The maximum 
reading using this probe was 18.7 kGauss = 1.87 Tesla (T). 
This value differs by less than one percent from the value 
stamped on the driver by the manufacturer: 1.858 T. We will 

therfore use a value of B=1.86 ± 0.01 T for this driver. We 
then counted the turns of aluminum wire on the voice coil (n 
= 30 ± l turns) and calculated 1 using the relationship: 

1 = n d n. ( I V-6 ) 
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For the 9.98 - 0.005 cm diameter coil we calculated a value 
of l=9.4±o.3m. Multiplying these two values together 
gives us for the transduction coefficient: 

B1 = 17.5 i 0.6 N/A. (IV-7) 

The other method used to determine B1 used the 
relationship : 

F = B1I = -kx or B1 = -kx/I , (IV-8) 

where I is the direct current applied to the coil, k is the 
spring constant measured earlier, and x is the displacement 
of the voice coil produced by the current. For this measure- 
ment on the 290-16K. driver we had used a Linear Variable 
Differential Transformer (LVDT). For the NIB driver, 
however, we used Mechanical Technology Incorporated’s MTI- 
1000 photonic sensor (serial 38010818) with plug - in 
cartridge model MTI-3812 (serial 1366) and probe model 125R. 
Figure IV-5 shows the photonic sensor with the NIB driver. 
One major advantage of the photonic sensor over the LVDT is 
that the photonic does not need to make physical contact 
with the surface it’s measuring (it uses light) whereas the 
LVDT does. With the Altec 290-16K’s rigid phenolic dia- 
phragm, the contact with the LVDT probably did not affect 
the measurement. The NIB diaphragm, however, is made of thin 
titanium and is quite flexible. It would have been impos- 
sible to get an accurate reading of the displacement of this 
diaphragm using the LVDT. We calibrated the photonic sensor 
using calibration fixture KD-CH-III A. The output voltage 
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Figure IV-5 . MTI-1000 Photonic Sensor with NIB Driver. 
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from the sensor was read on a HP 3478A multimeter (serial 
A110136). We then plotted DC voltage vs. displacement in 
millimeters. The resulting graph is shown in Figure IV-6 . 
The slope of this line gives us the calibration constant: 

dV/dx = -2.086 ± 0.001 V/mm. (IV-9) 

This value is within two percent of the calibration value 
provided by the manufacturer of dV/dx = 2.045 V/mm (back 
slope) . 

The next step was to run DC current through the coil 
and record the displacement of the diaphragm (via a voltage 
reading from the photonic sensor) . For the voltage supply we 
used the HP 467A power amplifier (serial 444-00243) . We 
found the current by measuring the voltage (with the HP 
3476A digital multimeter, serial 292162) across a ten ohm 

( 1 ^&) resistor placed in series with the circuit, and 
dividing the voltage by ten ohms. The resulting graph is 

shown in Figure IV-7. The slope of the line gives us : 

dl/dV = 2.10 ± 0.005 A/V. (IV-10) 

Combining the values in Equations IV-9 and IV-10 with the 
previously calculated high end value for the suspension 
stiffness k = 9.2 ± 0.6 X 10* kg/s 2 gj vcs : 

B1 = k (dx/dV) (dV/dl ) = 21.5 ± 1.4 N/A. (IV-11) 

This 15% difference between the values of B1 using 
the two independent methods is what prompted us to test for 
B1 and k using the static loading method. Our primary 
concern is that we would expect the direct measurements of B 
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CURRENT (A) DC VOLTAGE (V) 



PHOTONIC CALIBRATION 




Figure IV-6. Calibration Data for Photonic Sensor 
Voltage Output vs. Displacement. 



NIB - PHOTONIC CURRENT CALIBRATION 




VOLTAGE (V) 

Figure IV-7. Least Squares Fit of DC Current Through 
Voice Coil vs. Photonic Voltage Output. 
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and 1 to give us the highest possible value for B1 . This is 
because the maximum magnetic field may not be acting on the 
full length of the coil. We did not expect the value of B1 
calculated using the spring constant to be 15% higher than 
the direct measurement value. 

We were able to obtain static loading values for 
both B1 and k using the same equipment setup and the 
formulas : 

B1 = Amg/AI and k = - Amg/Ax (IV-12) 
For these measurements the voice coil is driven with a DC 
current. We then add previously measured masses (putty and 
soldering wire) to the voice coil, and measure the resulting 
displacement of the diaphragm using the photonic sensor. To 
measure B1 we determine how much current is needed to bring 
the voice coil back to its null displacement (original 
output voltage reading from the photonic sensor) after the 
mass is added. We measure k by keeping the current constant 
and measuring the difference in photonic output voltage 
(multiplied by the calibration constant to get displacement) 
resulting from the added mass. Table IV-1 lists the results 
of these measurements . 

The bracketed value in each k and B1 pair was 
calculated using the measured value in the pair and the 
relationship : 

BlI = -kx. ( IV- 13 ) 
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From these measurements we get the value k = 6.4 + o.2 kg/s = 
as reported earlier, and the transduction coefficient: 

B1 = 15.0 ± 0.2 N/A. ( IV- 14 ) 
For the model runs we used a value of B1 = 19.0 for the high 
k and B1 = 15.0 for the low k. 



TABLE IV-1. STATIC LOADING OF NIB DRIVER TO DETERMINE 
SUSPENSION STIFFNESS AND TRANSDUCTION COEFFICIENT 



Am ( g) 


AI (mA) 


AV (V) 


k (X10'* kg/s 2 ) 


B1 (N/A 


17 . 1 


-- 


5.2 


6 . 7 


[15.3] 


17.1 


-- 


5 . 5 


6 . 4 


[14.6] 


34 . 7 


o o ^ 


-- 


[6.4] 


14.6 


34.7 


22.3 




[6.7] 


15 . 3 


34.7 


— — 


. 288 


6 . 1 






(front slope) 



3 . Mechanical Resistance 

The next parameter we measured was the mechanical 

resistance (R m ) using the relationship: 

R m = 27tf 0 m 0 /Q. ( I V - 1 5 ) 

We attempted to find the driver quality factor (Q) using the 

same modified impedance technique described in Section 
III.G. As can be seen in Figure IV-8, the graph of the 
modified impedance magnitude vs. frequency is not a smooth 
resonance curve, as it was for the 290-16K (see Figure III- 
9). In fact, it looks like the superposition of two 
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NIB IMPEDANCE MAGNITUDE 




Figure IV-8. NIB Modified Impedance Magnitude 
vs. Frequency. 
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different impedance curves, the first one shorter and fatter 
and the second one taller and thinner. It was impossible to 
accurately measure half-power points from this curve, so we 
chose the phase angle method. (It is interesting to note 
that the JBL 2445J (ceramic) driver did have a smooth modi- 
fied impedance curve similar to that of the 290-16K.) Figure 
IV-9 shows the graph of phase angle vs. frequency over a 
wide range of frequencies. Although this curve generally 
approximates an arctangent, the segment of the line immedi- 
ately adjacent to the resonant frequency (assumed to be the 
frequency where 0 = 0 , in this case 534.5 ± 0.5 Hz) is linear 
(for a frequency range of about ten Hz) . This segment is 
plotted in Figure IV-10. The slope of the line is de/df = - 

42.0 ± 0.5 X 10 _3 s (after converting degrees to radians) . As 

in Section II.G, Q is given by: 

0 = -(f„/2) (de/df ) = 11.2 ± 0.2. ( IV- 16 ) 

This gives a value for R m (using the moving mass in air 
m <=. = 7.3 ± 0.1 g ) : 

R m = 2.34 ± 0.05 kg/s. ( IV- 1 7 ) 

We decided that a more accurate measurement of Q and 

R m could be found by making the Q measurement in vacuum. For 
this we used the open circuit free decay method. This method 
consists of sending a signal to the driver at resonance with 
an HP 3314A function generator (serial A342044) , then 
simultaneously switching off the signal and triggering a 
Nicolet model 310 storage oscilloscope (serial 88004641) to 
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PHASE ANGLE (DEG) PHASE (DEG) 




Figure IV-9. NIB Modified Impedance Phase Angle 
as a Function of Frequency--Full 
Frequency Range. 



NIB PHASE ANGLE - LINEAR SEGMENT 




Figure IV- 1 0 . NIB Modified Impedance Phase Angle 

as a Function of Frequency--Linear 
Segment Near Resonance. 
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record the free decay. A sample of the output from this 
process is shown in Figure IV-11. The equation for a simple 
harmonic oscillator experiencing free decay is: 

A(t) = A 0 exp (-t/t) sin(ie«t+#) , (IV-13) 

where A(t) is the amplitude as a function of time, A 0 is the 
initial (maximum) amplitude, t is the time (= nT , or the 
number of periods times the period) , x is the free decay 
constant (= QT/jt, see Equation III-29) , oj 0 is the angular 
resonant frequency, and is the phase angle. We recorded 
the peak amplitude as a function of period number, and 

plotted the natural logarithm of the peak amplitude vs. 
period number. The equation representing this line is: 

In A ( T ) = In Ao - nT/ t- . (IV-19) 

The slope of this line is t/x = nT/t = n?t/Q. Therefore Q can 

be found by dividing ft by the slope of the line. Figure IV- 

12 shows a semi-log plot for the NIB driver in air. The 

value of Q for this plot is 21.5, which is almost twice the 

value we found using the modified impedance phase angle 

measurements (Q = 11.2). Obviously the air strongly affects 
the measurements. The scattering of the data points in this 
figure is probably due to the same phenomena that causes the 

double peaks in the impedance curve. Figure IV-13 shows the 

semi-log plot of the amplitude vs. period number for vacuum. 
This curve is much smoother and, we assume, more reliable. 
The Q for this data is: 

Q = TT/slope = 19.8 ± 0.05. (IV-20) 
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Figure IV-ll. Sample of Free Decay Output from 

Nicolet 310 Storage Oscilloscope. 
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In amp 



NIB FREE DECAY - AIR 




Figure IV-12. Semi-log Plot of Amplitude 

vs. Period Number for NIB 
Free Decay in Air. 



NIB FREE DECAY -VACUUM 




Figure IV-13. Semi-log Plot of Amplitude 

vs. Period Number for NIB 
Free Decay in Vacuum. 
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We can also find the period from the free decay curve (and 
thus the resonant frequency as the inverse of the period) . 
We use the oscilloscope to measure the time for ten cycles 
and divide this number by ten. For the vacuum run the 
resonant frequency was 839 ± 0.5 Hz. These values (and the 
moving mass = 3.2 g in vacuum) give us the mechanical resis- 
tance in vacuum: 

Rm = 27T-f 0 nWQ = 0.83 ± 0.01 kg/s. (IV-21) 
This is the value we used for the model. 

Table IV-2 gives a comparison of applicable 
parameters for the three drivers. Notice that the parameter 
Bl/R m * = ^,m 0 is almost three times greater for the NIB driver 
than the 290-16K. This is the parameter that we said (Sub- 
section IV. A. 2) we wanted to maximize to get the highest 
thermoacoustic efficiency. Thus the driver that gives us the 
best height and weight characteristics (see photograph of 
the three drivers in Figure IV-14) also gives us the best 
electroacoustic efficiency. 

C. MODELED EFFICIENCY 

In order to model the electroacoustic efficiency of the 
driver, we need to know what acoustic load (impedance) it 
will see. The final resonator and stack composition will not 
have been determined by the time this thesis is completed. 
We therefore assumed a range of possible values for the 
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TABLE I V-2 . COMPARISON OF PARAMETERS FOR 
THE THREE DRIVERS 



Property (units) 


Altec 290 


JBL 2445 J 


JBL NIB 


Moving mass, m 0 (kg) 


0 . 0151 


0 .015 


0.015 


Spring constant (N/m) 


3 . 52X1 O'* 


~8 . 5X10** 


-9.0X10 4 


Magnetic Induction, 
B (Tesla) 


2 . 1 


1.97 


1 . 87 


Voice Coil 

Diameter (cm) 


7.22 


9.98 


9 . 98 


B1 (N/A) 


21 . 5 


18.5 


17.5 


Voice coil DC 
resistance ( ?i ) 


10 . 6 


7.1 


7.6 


Mechanical 

resistance (kg/s) 


2 . 77 


**1 . 0 


0 . 83 


Voice coil 

inductance (H) 
(at 20 kHz) 


1 . oxio-* 


1 . 14X1 0 - "* 


1 . 14X10- 


Mass (kg) 


14.5 


13.8 


3 . 1 


Height (in) 


6.0 


5.5 


2 . 5 


Magnetic energy 
density (J/cm 3 ) 


~0 .5 


■vO . 5 


3.2 


Bl/R^^c^nic 


510 


1230 


1410 



( m/Coul-kg ) 
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Figure IV-14. Photographic Comparison of the 

Three Drivers — Altec 290-16K, 
JBL 2445J , and JBL NIB. 
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acoustic load based on the experimental data acquired to 
date. 

The model parameters that relate to the acoustic 
load are the density of the gas (RHO) , the speed of sound in 
the gas (C), the resonator quality factor (Q) , the resonator 
resonant frequency (FO), and the resonator length (L). The 
model was designed for a closed-ended tude of uniform cross- 
section. The STAR resonator has an open tube (with a bulb 
attached) and nonuniform cross-section. The problem of 
applying the existing computer model to the STAR resonator 
was solved by scaling the STAR parameters to make it look 
like a closed-ended uniform cross-section tube. 

1 . Input Values for Pure Helium 

For pure helium, the experimental data available was 
the range of resonator quality factor values (10-20), and 
the speed of sound and gas density at room temperature 
(about 300K) and ten atmospheres (c RT = 1024 m/s, 

= 1.618 kg/m 3 ). Also available were four resonant fre- 
quencies, two for a short resonator at warm (580 Hz) and 

cold (520 Hz) temperatures, and two for a long resonator at 
warm (520 Hz) and cold (470 Hz) temperatures. Using this 
information, we first calculated the effective lengths of 
the resonators at room temperature as if they were closed- 
ended, uniform cross-section tubes. For this we used the 
relationship : 

A = 2L = c/f, or L a<lf = c/2f , (IV-22) 
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where h is the wavelength, L is the length of the tube, c is 
the speed of sound, and f is the resonant frequency. This 
calculation gave us effective lengths for the resonators: 
L* = 0.883m and Li = 0.985m. For the computer runs we use a 
convenient value of = 0.90m. We then calculated the 
speed of sound for the operating temperature of approxi- 
mately 250K using the relationship: 

c ( T ) = c R t [T/Tr T ] 1/2 . (IV-23) 
This gives us a speed of sound c =!= o = 935 m/s. If we use 
this sound speed and the effective resonator lengths to 
estimate the resonant frequencies for 250K , we get f m = 529 
Hz and f x = 474 Hz, which agree with the measured values. We 
calculated the gas density at the operating temperature 
using the relationship: 

p = tfP/c 2 = 1.94 kg/m 3 , ( IV-24 ) 
where % is the ratio of the specific heats (1.667 for 
helium) and P is the pressure. 

Table IV-3 lists the input parameters for the pure 
helium computer runs. Note that for the computer model, the 
moving mass must take the mass of the reducer cone and 
bellows into account (see Figure IV-15) . The total moving 
mass is expected to be approximately 15 g. 

2 . Input Values for Helium-Xenon 

There was somewhat less data available for the 
helium-xenon mixture (12.57° Xe) . We had the sound speed at 
300K and ten atmospheres (c R t = 472 m/s). We also had 
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TABLE I V-3 . COMPUTER MODEL INPUT-- 
NIB DRIVER WITH HELIUM AT TEN ATMOSPHERES 
AND TEMPERATURE = 250K 



Symbol 


Meaning 


Units 


NIB 


RM 


Mechanical resistance 


kg/s 


0 . 83 


M 


Moving mass 


kg 


0 .015 


SSUP 


Spring constant 


N/m 


6.6, 9.0 X10 


RHO 


Density of gas 


kg/m 3 


1 . 94 


C 


Sound speed 


m/s 


935 . 0 


S 


Resonator X-sec . area 


m 2 


1 . 14X1 0 -3 


BL 


Driver transduction factor 


N/A 


15.0, 19. 


Q 


Resonator quality factor 




10. 15,20 


FO 


Resonator resonance freq. 


Hz 


519.4 


A 


Back volume X-sec. area 


m 2 


4. 42X10" 3 


VOL 


Back volume 


m 3 


3. 28X10-* 


LE 


Voice coil inductance 


H 


1 . 14X10-* 


RE 


Voice coil DC resistance 


a 


7 . 6 


I 


Driver current 


A 


1 . 0 


L 


Resonator length 


m 


0.90 
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(a) 




i 

i 

i 



I 

i 



(b) 

Figure IV-15. Photograph of JBL NIB Driver 

with (a) Reducer and (b) Bellows. 
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resonant frequencies for the same four conditions as for the 
pure helium: short-warm (265 Hz) , short-cold (230 Hz) , long- 
warm (235 Hz), and long-cold (215 Hz). We used Equation 
IV-22 to find the effective resonator lengths (L„ = 0.391m 
and L i = i.004m). However, once again we used a convenient 
effective length of 0.90 m for the computer runs. We used 
Equation IV-23 to find the sound speed for He-Xe at the 
operating temperature, = 431 m/s. The estimated reso- 
nant frequencies for the cold temperature were f m =242 Hz and 
fi=2l5 Hz, once again comparing favorably with the measured 
values. To find the density of the mixture at the operating 
temperature we used a form of the ideal gas law: 

p = PM/RT, ( IV-25 ) 
where M is the atomic mass of the mixture, R is the 
universal gas constant, and P and T are the pressure and 
temperature. We calculated the atomic mass of the mixture 
using the equation: 

Mm = Mh«Xh. + Mx-Xx. . (IV-26) 

where the M's are atomic masses (4 for helium and 131.3 for 

xenon) and the x's are percentages of the gases in the 
mixtures. This equation gives us an atomic mass for the 
12.5% Xe mixture of 19.9 g/mol , and a density of 9.70 kg/m 3 
at 2S0K and ten atmospheres. Table IV-4 lists the computer 
input values for He-Xe, including a range of Q's from 5 to 
30. (We are not sure what to expect for the actual value of 
Q. since a new stack will also be used.) 
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TABLE I V - 4 . COMPUTER MODEL INPUT-- 
NIB DRIVER WITH HELIUM-XENON AT TEN ATMOSPHERES 

AND TEMPERATURE = 250K 



Symbol 

RM 

M 

SSUP 

RHO 

C 

s 

EL 

Q 

FO 

A 

VOL 

LE 

RE 

I 

L 



Meaning 

Mechanical resistance 

Moving mass 

Spring constant 

Density of gas 

Sound speed 

Reso: .or X-sec . area 

Driver trans. factor 

Resonator quality factor 

Resonator resonance freq. 

Back volume X-sec . area 

Back volume 

Voice coil inductance 

Voice coil DC resistance 

Driver current 

Resonator length 



Units 

kg/s 

kg 

N/m 

kg/m 3 

m/s 

m = 



NIB 
0 . 83 
0 . 015 

6.6. 9.0 X10- 
9 . 70 
431 . 0 
1 . 14X10“° 



N/A 15.0. 19.0 

5,10,15,20,25,30 



Hz 



m 3 



m 3 

H 



A 

m 



239 . 4 
4.42X10-° 
3 . 28X1 0 - "* 
1 . 14X10—* 
7.6 
1.0 
0.90 
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3 . Computer Model Results 



Figures IV-16 through IV-19 are sample output plots 
for pure helium with B1 = 15 N/A, k = 6.4 X10-* kg/s 2 , and 
Q = 20. Notice that the peak values for the driver and 
resonator are almost matched (the driver resonance is at the 
lower frequency) . The thermoacoustic efficiency for this set 
of values is about 70.4%. Figures IV-20 through IV-23 show 
the corresponding plots for the helium-xenon mixture. This 
time the peaks are not matched (the driver peak is off the 
scale) and the electroacoustic efficiency is about 16.6%. 
The primary reason for the difference between the two gases 
is the lower sound speed, and therefore lower resonant fre- 
quency, in the He-Xe mixture. 

Table IV-5 lists the maximum value of 
electroacoustic efficiency for each of the test runs. The 
consistently low efficiencies for the He-Xe emphasize the 
importance of matching the two resonant peaks. We recommend 
that the driver frequency be adjusted to match the resonator 
operating frequency. This can best be done by reducing the 
stiffness of the titanium surround (see Section III.H). 
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Figure IV-16. Computer Model Output for NIB with 

Helium ( B1 = 15, Q = 20)--Driver 
Voltage vs. Frequency for One Amp. 
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Figure IV-17. Computer Model Output for NIB with 

Helium (Bl = 15, O = 20)--Piston 
Velocity vs. Frequency for One Amp. 
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Figure IV-18. Computer Model Output for NIB with 

Helium (B1 = 15, Q = 20 ) --Acoustic 
Power vs. Frequency for One Amp. 
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Figure IV-19. Computer Model Output for NIB with Helium 

(B1 = 15, Q = 20 ) --Electroacoustic 
Efficiency vs. Frequency. 
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Figure IV-20. Computer Model Output for NIB with Helium- 

Xenon (B1 = 15, Q = 20)--Driver Voltage 
vs. Frequency for One Amp. 
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Figure IV-21 . Computer Model Output for NIB with Helium- 

Xenon (B1 = 15, Q = 20)--Piston Velocity 
vs. Frequency for One Amp. 
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Figure IV-22. Computer Model Output for NIB with Helium- 

Xenon (B1 = 15, Q = 20 ) --Acoustic Power 
vs. Frequency for One Amp. 
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Figure IV-23. Computer Model Output for NIB with Helium- 

Xenon (B1 = 15, Q = 20 ) --Electroacoustic 
Efficiency vs. Frequency. 
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TABLE IV- 5 . COMPUTER MODEL RESULTS FOR NIB DRIVER- - 
PEAK ELECTROACOUSTIC EFFICIENCY (%) 



Q 


He ( Bl=15 ) 


He ( Bl=19 ) 


He-Xe ( Bl=15 ) 


He-Xe ( Bl = 19 


5 






5.1 


5 . 9 


10 


71 . 3 


70 . 4 


9.2 


10 . 7 


15 


71.6 


74.3 


13. 1 


15.1 


20 


70 . 4 


76.2 


16.6 


19.1 


25 






19.6 


22.7 


30 






22 . 8 


26 . 0 
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V. STAR DRIVER CONSTRUCTION 



A. INTRODUCTION 

This chapter describes the construction of the STAR 
driver. A large part of this thesis involved the design of 
the parts necessary to integrate the NIB driver with the 
resonator and GAS cannister, as illustrated in Figure 1-4. 

The major design work was for the driver housing. The 
requirements for the driver housing are: 

1 - lightweight, 

2 - strong (pressurized to 210 psi), 

3 - interface to GAS canister lid (12 10-32 bolt holes 

centered at 9.5 inches diameter), 

4 - interface with resonator (preferably interchangeable 

with earlier test resonators) , 

5 - helium leak-tight feed-through ports for electrical 

leads (driver electric power, accelerometer and 

microphone outputs) , and 

6 - hold NIB magnet structure securely. 

Figure V-l shows a cross-sectional drawing of the driver 
housing (not including any of the feed-through ports, which 
will be discussed later). Distinction between parts is shown 
by different hatching. Part A is the NIB magnet structure 
which was discussed in Chapter IV. Part B is the main body 
of the driver housing, machined out of 6061-T6 aluminum. 
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Figure V-l. Cross-sectional Schematic of the STAR Driver 
Housing Showing (A) NIB Driver, (B) Driver 
Housing Main Body, (C) Aluminum Pusher Plate, 
(D) Pressure Lid, (E) Leak Test Lid, 

(F) Resonator Vacuum Housing, and (G) Bolt 
Holes for Interface with Resonator. 
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Part C is a 6061-T6 aluminum pusher plate used to hold the 
magnet structure securely in the body using six 10/32 bolts. 
Part D is the pressure lid, sealed a l/16th inch diameter 
lead o-ring and held in place with six 10/32 bolts. Part E 
represents part of the leak test lid (see Figure C-l). When 
completed this will be replaced by the GAS can lid provided 
by NASA. The driver housing must match up with the GAS can 
lid’s 9.5 inch diameter 10-32 bolt circle. Part F is the 
resonator section housing. The arrows labeled G point to the 
bolt holes which form the driver-resonator interface. 

B. DESIGN CONSIDERATIONS 

The driver-resonator system will be pressurized to ten 
atmospheres with a helium-xenon mixture (12.5% xenon). 
NASA’s regulations say that it must be pressure tested to 
1.5 times the operating pressure (NASA, 1986a), or 15 atmo- 
spheres (=1.52X10 fe pa = 220 psi). An added complication of 
working with helium (or helium mixtures) is that it diffuses 
through almost all materials except metals. As a result, we 
are not able to use rubber o-rings to seal metal parts 
together, but must use lead o-rings instead. 

1 . Electrical Feed-Throughs 

NASA’s regulations prohibit us from leaking gases 
into space, primarily to avoid contamination of other exper- 
iments and for the safety of the Shuttle and astronauts. 
This regulation is important not only for the major o-ring 
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seals , but also for the microphone and accelerometer outputs 
to the frequency control electronics (see Figure 1-5), and 
connecting the audio amplifier to the driver. There are 
three individual wires and three co-axial connectors that 
must penetrate the driver housing without allowing leakage. 
The other housing penetrations are a DC pressure gage for 
the driver housing and a fill line from a helium-xenon 
bottle which will stabilize the pressure in the system 
against temperature changes in the environment and those 
generated by the refrigerator. 

Based on prior experience in other applications, we 
decided that Emerson & Cuming’s STYCAST <tm> 2850FT epoxy 
(Canton, MA 02021) would be the best choice for providing 
the electrically insulating seal for the wires in the 
housing. STYCAST <tm> 2850FT is a quartz-filled epoxy formu- 
lation with remarkably good over-all general properties. In 
addition to having excellent electrical grade insulation 
properties, good adhesion, and excellent resistance to chem- 
icals and solvents, 2850FT has unusually high thermal con- 
ductivity and low thermal expansion. This combination of 
properties makes this epoxy valuable for applications where 
electrical insulation and mechanical protection must be 
maintained while coping with heat transfer considerations. 

One definition for a sealed container is a helium 
leak rate of less than 10~ 4 standard cc ’ s per second (NASA. 
1986b) . To test the epoxy we designed and built a test 
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fixture to be used with the Veeco MS-9 helium leak detector. 
A cross-sectional diagram of the test fixture is shown in 
Figure V-2. Part A is a copper tube with 1-1/8" outer diam- 
eter to interface with the leak detector. Soldered to the 
copper tube is a piece of solid brass stock about one inch 
thick (B) with a 0.052 inch diameter hole bored through the 
center. Epoxied into the hole is a 0.025 inch diameter bare 
copper wire (C) . This design is motivated by the desire to 
provide as long a travel path and as small a cross-section 
as possible for the helium gas molecules. Part D is a hollow 
aluminum housing designed to hold pressurized helium. We 
decided that a rubber o-ring should hold the helium suffi- 
ciently well for a short term leak test. Part E is a 1/8” 
NPT (National Pipe Thread) hole for connecting to a helium 
fill line. A photograph of the test fixture with and without 
the housing is shown in Figure V-3. 

The copper tube and brass stock were soldered toge- 
ther. We boiled the assembly to remove excess particulates 
and rinsed it with alcohol. We then etched the hole in the 
brass and the copper wire with nitric acid diluted with 
distilled water (about 20% acid) to provide a clean, rough 
surface for the epoxy. Before etching, cellophane tape was 
placed on the portion of the brass stock where the o-ring 
would be sealing (after making sure that the tape would not 
dissolve in the acid) . We then rinsed the assembly with 
distilled water and dried it, mixed the epoxy, and inserted 
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Figure V-2. Cross-sectional Schematic of the Helium 
Leak Detector Test Fixture Showing 

(A) Copper Tube for Leak Detector Interface, 

(B) Solid Brass Stock, (C) Bare Copper Wire, 
(D) Aluminum Housing, and (E) Connection for 
Helium Fill Line. 
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Figure V-3. Photograph of the Helium Leak Detector 
Test Fixture. 
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the wire into the hole in the brass stock, making sure the 
epoxy filled the gap. 

The following procedure was used to run the leak 
detector tests: first zero the meter, then do a test run on 
the sensitivity calibrator. The sensitivity calibrator used 
was a Veeco type SC-4 serial number 3179 with a leak rate of 

6.5X10 _e std cc/sec (standard cubic centimeters per second). 
We would then remove the calibrator, recheck the meter for 
zero, attach the test fixture, read the meter, remove the 
fixture, and recheck for zero. 

The test fixture was initially pressurized to 115 
psi with helium. We tried to test it immediately but found 
that the detector would not draw enough of a vacuum to check 
for helium. We thought that there might be outgassing from 
flux or something in the tube. So we pulled a vacuum on the 
tube overnight and the next day it worked fine. We tested 
the fixture on the first, third, and sixth days after it had 
been pressurized. The numerical results are listed in Table 
V— 1 . Even after we had waited six days for the helium to 
diffuse through the epoxy, the fixture leakage rate was an 
order of magnitude less than the calibrator. The leakage 
rate was therefore four orders of magnitude less than that 
allowed by NASA to be considered leak-tight (1X10-^- s td 
cc/sec). After the third leak test we checked the pressure 
in the aluminum chamber and found it to be within a few 
percent of the original pressure. We therefore concluded 
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that the bare wire-epoxy method with a length/diameter ratio 
of about 20 could be used for the driver housing penetrating 

electrical feed-throughs. 

TABLE V-l. HELIUM LEAK DETECTOR TEST RESULTS 
ON COPPER WIRE-EPOXY FIXTURE. 



Number days after 


null meter 


calibrator 


copper wire feed 


pressurization 


reading 


reading 


-through reading 


1 


0 


1 . 1 


0 . 02 




0 


1 . 2 


0 . 08 


6 


0 


1 . 1 


0.10 



Figure V-4 shows a close-up of the feed-through plug 
for the wire leads. It will house a total of five wires 
(with space available for a sixth if necessary) . Two of the 
leads will be for the voltage supply to the voice coil. 
These leads will consist of bare copper wires epoxied into 
the small holes seen in the figure. Each wire will then be 
soldered to a 6/32 phillip-head screw with a 1-72 blind- 
tapped hole in it. The screw will be glued into an 8/32 
hole, forming a mechanically stable, demountable, electrical 
terminal. The accelerometer and microphone leads will also 
have bare copper wire epoxied into the small hole, and will 
then be connected to Microdot ‘™> coaxial fittings (Malco 
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Figure V-4 . Driver Housing Feed-Through Plug 
for the Electrical Leads. 
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part number 031-0001—0001, South Passadena, CA 91030} . 
Figure V-5 shows an example of each of these leads . 

2 . Measurement Devices 

The microphone and accelerometer used by STAR are 
similar to those used by Hofler (1986 and 1988). In fact, 
the resonator end of the driver housing was designed to 
accept Hofler's resonator, calibration volume, and instru- 
mentation for testing. The microphone used to measure the 
acoustic pressure in the resonator is a small quartz-crystal 
( Valpey-Fisher , Hopkinton, MA 01748) or PZT dynamic pressure 
transducer (Piezoelectric Products lot number 4-9-84, 
F013184, Metuchen, NJ). It is mounted in the driver housing 
where the housing connects to the resonator flange. The 
microphone will be calibrated by comparison to a standard 
transducer, which will have been previously calibrated via a 
reciprocity technique. The microphone will be connected to a 
miniature FET follower circuit. The FET buffer (Eletec model 
304 impedance converter, Daytona Beach, FL 32020) converts 



the microphone's high impedance 


output signal 


to a 


low 


impedance 


to 


reduce 


sensitivity 


to electrical 


noise 


and 


maintain 


the 


output 


voltage and 


phase while 


driving 


the 



cables necessary to reach the signal processing electronics. 
The buffer's DC power supply and output cables are complete- 
ly shielded through the driver housing to the buffer (see 
Figures C-3 and C-7 for location of microphone and FET in 
driver housing) . 



139 




(a) 




(b) 

Figure V-5. Example of Wire Lead Connectors 

for Feed-Through Plug: (a) Single 

Wire and (b) Co-axial Cable. 
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The Endevco model 22 Picomin < " r ' v,> accelerometer is 
attached to the side of the bellows stiffener facing the 

driver. The integrated accelerometer signal will be used as 
a measure of the bellows velocity, which will be calibrated 
by driving a small cavity of calculable acoustic impedance. 
Thus the product of the outputs from the microphone and 
integrated accelerometer can be combined to yield the acous- 
tic power and their relative phase can be used to keep the 
driver frequency at the tube resonance. The operation of the 
refrigerator at resonance is essential both to optimize the 
acoustic power delivered to the "stack” and to guarantee 
that the "stack” is always located at the proper position in 
the standing wave as discussed in Sub-section II. A. 4. 

The DC pressure transducer is an Omega <TM> PX-80 
designed for high pressure applications requiring long-term 
stability. Figure V-6 shows a drawing of the pressure trans- 
ducer and its driver housing fitting. 

The helium port in the housing is provided to 
empty and fill the driver and resonator and so that a stable 
pressure may be maintained in the driver and resonator 
through the ballast volume provided by the reservoir (see 
Figure 1-4 ) . 

There is a 1/32 inch diameter hole filled with a ten 
thousandths inch inner diameter half inch long capillary in 
the driver housing where it connects to the resonator. The 
function of the capillary is to equilibriate the static 
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Figure V-6. 



DC Pressure Transducer and Driver 
Housing Plug Fitting. 



pressure between the driver and resonator during filling and 



non-operating times while maintaining a pressure seal 
acoustic operating frequencies. 

A complete set of dimensioned drawings for 
driver housing and its auxiliary equipment is given 
Appendix C. A list of materials for the driver portion 
the STAR is given in Appendix D. 



at 

the 

in 

of 
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VI . CONCLUSION AND RECOMMENDATIONS 



The objective of the STAR project is to test and space 
qualify a new continuous cycle cryogenic refrigeration 
system for cooling of sensors and electronics which is based 
upon the newly-discovered thermoacoustic heat pumping 
effect. The new refrigerator has no sliding seals, a cycle 
frequency of about 300 Hz, and uses acoustic resonance to 
enhance overall efficiency. This thesis is concerned speci- 
fically with the design and testing of the electrodynamic 
transducer which is responsible for the electro-acoustic 
power conversion. A computer model of the driver/resonator 
system was presented along with the techniques for measure- 
ment of the electrical and mechanical parameters used as 
input for the model . 

A final driver design (including dimensional drawings) 
utilizing a modified JBL 2450J neodymium-iron-boron 
compression driver and associated leak-tight electrical 
feed-throughs, microphone, accelerometer, pressure gage, 
pressure housing, and resonator interface is provided. 

Since we only received the NIB driver in late December 
1987, there wasn’t time to perform enough parameter 
measurements to obtain accurate values for the suspension 
stiffness (k) and transduction coefficient ( B1 ) . The next 
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group of students should find these values for use in the 
model and efficiency calculations. 

The resonator experiments (not reported here) indicate 
that Helium-Xenon may be the preferred gas mixture for the 
STAR. If this is the case, the NIB driver resonant frequency 
will have to be lowered for it to match the operating 
(resonant) frequency of the resonator. This could be 
accomplished by adding moving mass to the voice coil- 
reducer-bellows assembly, but that would decrease the 
velocity and lower the electroacoustic efficiency. The 
preferred option for driver resonance adjustment is to 
decrease the suspension stiffness by cutting away part of 
the surround (see Section III.H). Since the NIB voice coil 
has a titanium surround and diaphragm this operation will 
have to be performed carefully. The titanium diaphragm 
(dome) will also have to be removed. (The dome's function is 
performed by the reducer and bellows.) 

The task for the next group of students is to assemble 
the NIB driver with its housing and accessories (feed- 
throughs, microphone, accelerometer, etc.) and perform 
efficiency measurements on the ground similar to those that 
will be performed in space. This series of measurements will 
closely resemble those performed on the Altec 290- 16K driver 
as described in Chapter III. 
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APPENDIX A. COMPUTER MODEL CALCULATION PROGRAM 



PROGRAM FRIDGE 

C INTRODUCTION: 

C 

C THE PURPOSE OF THIS PROGRAM WRITTEN BY LT . MICHAEL P. 

C SUSALLA IS TO CALCULATE THE DATA POINTS AND PLOT THE 
C PERFORMANCE CHARACTERISTICS OF A CONCEPTUAL MODEL OF A 
C THERMOACOUSTIC HEAT ENGINE. THE COMMENTS INCORPORATED INTO 
C THIS PROGRAM AND THE COMMENTS INCORPORATED IN THE PLOTTING 
C PROGRAM "GRAFRI FORTRAN Al", WILL ALLOW THE USER TO PRO- 
C DUCE FOUR GRAPHICAL PLOTS WHICH WILL CLOSELY APPROXIMATE 
C THE PERFORMANCE OF AN ACTUAL THERMOACOUST IC HEAT ENGINE. 

C THE GOAL OF THE PROGRAM /SIMULATION' 'S TO ALLOW THE 

C OPERATOR TO VARY DIFFERENT PARAMETERS TO ACHIEVE WHAT IS 
C BELIEVED TO BE AN OPTIMUM CONDITION WHERE THE ENGINES 
C DRIVER RESONANCE IS MATCHED WITH THE TUBE RESONANCE. 

Q 

C GENERAL INSTRUCTIONS: 

C 

C THIS PROGRAM IS SET UP TO PROVIDE FOUR PLOTS: VOLTAGE 

C VERSUS FREQUENCY, PISTON VELOCITY VERSUS FREQUENCY, ACOUS- 
C TIC POWER VERSUS FREQUENCY, AND ELECTROACOUSTIC EFFICIENCY 
C VERSUS FREQUENCY. IN GENERAL THE PROGRAM HAS HARD INPUT 
C DATA SUCH THAT ALL YOU MUST DO TO GET A SET OF PLOTS IS 
C COMPILE "GRAFRI FORTRAN Al" IN FORTVS AND THEN RUN "FRIDGE 
C FORTRAN Al" IN WF77 FROM FLIST ON THE NPS MAINFRAME USING 
C A TEKTRONICS 618 TERMINAL. THIS ACTION WILL GIVE YOU A SET 
C OF PLOTS FOR THE INPUT DATA. TO TRY OTHER DATA YOU MUST GO 
C TO XEDIT ON THE "FRIDGE FORTRAN Al" PROGRAM AND VARY THE 
C INPUT DATA. AS LONG AS YOUR CHANGES ARE RESTRICTED TO THE 
C "FRIDGE” PROGRAM THE ONLY ACTION NECESSARY TO PRODUCE NEW 
C PLOTS IS TO RUN THE PROGRAM IN WF77. IF YOU FIND IT NECES- 
C SARY TO CHANGE THE PLOTTING PROGRAM, "GRAFRI", FOR WHAT- 
C EVER REASON, IT MUST BE RECOMPILED IN FORTVS PRIOR TO THE 
C NEXT ATTEMPT TO PLOT OR THE CHANGES WILL NOT BE SEEN BY 
C THE MACHINE. 

C 

C INPUT SECTION: 

C A. DEFINITION OF VARIABLES: 

C ALL OF THE VARIABLES HAVE BEEN ESTABLISHED AS REAL 

C OR COMPLEX QUANTITIES BASED ON THE REQUIREMENTS OF THE 
C ARITHMETIC ALGORITHOM OF THE CONCEPTUAL MODEL. EACH OF 

C THE VARIABLES WILL NOW BE DEFINED. 

C LCOI L=THE IMAGINARY PART OF THE IMPEDANCE OF THE VOICE 

C COIL OF THE DRIVER, AND IS CALCULATED TO BE W*LE. 

C K=THE PROPAGATION CONSTANT OF THE WAVES IN THE TUBE 

C AND IS GIVEN BY W/C . 
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C IMZM=THE IMAGINARY PART OF THE MECHANICAL IMPEDANCE OF 

C THE TUBE/DRIVER COMBINATION. IT IS GIVEN BY A RATHER LONG 
C AND COMPLEX EQUATION IN THE CALCULATION SECTION OF THE 
C PROGRAM. 

C MAGU=THE MAGNITUDE OF THE COMPLEX QUANTITY "U" . IN 

C METERS/SEC. 

C IMZT=THE IMAGINARY PART OF THE TUBE MECHANICAL IMPE- 

C DANCE AND IS GIVEN BY SUBTRACTING THE DRIVER PORTIONS FROM 
C IMZM . 

C M=THE MOVING MASS OF THE SYSTEM, AND IS THE SUM OF THE 

C VOICE COIL AND THE DIAPHRAGM MASSES OF THE DRIVER, IN 
C KILOGRAMS. 

C L=THE LENGTH OF THE TUBE IN METERS. 

C LE=THE INDUCTANCE OF THE VOICE COIL, IN HENRYS. 

C I =THE INPUT CURRENT FROM THE SOURCE, IN AMPS. 

C MAGV=THE MAGNITUDE OF THE COMPLEX QUANTITY "V" , IN 

C VOLTS. 

C ZM=THE MECHANICAL IMPEDANCE OF THE TUBE/DRIVER SYSTEM, 

C AND IT IS MADE UP OF THE REAL PART, REZM , AND THE I MAGI - 
C NARY PART, IMZM. 

C ZE=THE ELECTRICAL EQUIVALENT OF THE MECHANICAL IMPE- 

C DANCE OF THE TUBE/DRIVER SYSTEM. IT IS GIVEN BY TRANSFOR- 
C MING ZM ACROSS A TRANSFORMER USING THE BL COEFFICIENT OF 
C THE DRIVER. 

C ZCOIL=THE COMPLEX ELECTRICAL IMPEDANCE OF THE VOICE 

C COIL. 

C ZTOT=THE TOTAL ELECTRICAL IMPEDANCE OF THE SYSTEM, 

C COMPLEX. 

C V=THE COMPLEX VOLTAGE OF THE SOURCE NECESSSARY TO 

C PRODUCE THE SELECTED INPUT CURRENT. 

C U=THE VELOCITY OF THE PI STON ,( DIAPHRAGM ) , OF THE 

C DRIVER IN M/S. 

C ZTUBE=THE COMPLEX MECHANICAL IMPEDANCE OF THE TUBE 

C ITSELF. 

C P=THE COMPLEX PRESSURE OF THE SYSTEM, IN NEWTONS/M ** 2 . 

C PI =A CONSTANT FOR CONVERSIONS TO AND FROM ANGULAR 

C UNITS. 

C RM=THE MECHANICAL RESISTANCE OF THE MOVING MASS OF THE 

C SYSTEM. 

C SSUP=THE MECHANICAL STIFFNESS OF THE SYSTEM DUE TO THE 

C DRIVER SUSPENSION. 

C RHO=THE DENSITY OF THE GAS IN THE TUBE AT A SPECIFIED 

C MEAN TEMPERATURE. 

C C=THE SPEED OF SOUND FOR THE SPECIFIED GAS AND 

C TEMPERATURE. 

C S=THE CROSS SECTIONAL AREA OF THE TUBE. 

C BL= A COEFFICIENT OF THE DRIVER THAT IS BASED ON THE 

C STEADY AIR GAP FLUX DENSITY OF THE DRIVER AND THE LENGTH 
C OF ITS VOICE COIL. 

C Q=THE QUALITY FACTOR OF THE TUBE. 

C FO=THE RESONANCE FREQUENCY OF THE TUBE. 

C A=THE CROSS SECTIONAL AREA OF THE BACK VOLUME OF THE 
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C DRIVER. 

C VOL=THE VOLUME OF THE BACK VOLUME OF THE DRIVER. 

C RG=THE INTERNAL ELECTRICAL RESISTANCE OF THE POWER 

C SOURCE. 

C RE=THE ELECTRICAL RESISTANCE OF THE VOICE COIL. 

C SVOL=THE MECHANICAL STIFFNESS OF THE SYSTEM DUE THE 

C BACK VOLUME OF THE DRIVER. 

C WO=THE RESONANCE FREQUENCY OF THE TUBE IN ANGULAR 

C UNITS, RAD/SEC. 

C ALPHA=A CHARACTERISTIC QUANTITY OF THE TUBE THAT IS A 

C FUNCTION OF THE RESONANT FREQUENCY. THE Q AND THE SPEED OF 
C SOUND. 

C F=FREQUENCY . IN HERTZ. 

C W=THE FREQUENCY CONVERTED TO RAD/SEC. 

C FORCE=THE FORCE ON THE MOVING MASS DUE TO THE INTER- 

C ACTION OF THE VOICE COIL MAGNETIC FIELD AND THE DRIVER 
C MAGNET. 

C REZT=THE REAL PART Of THE COMPLEX MECHANICAL TUBE 

C IMPEDANCE. 

C PU=THE PRESSURE TIMES THE VOLUME VELOCITY, ACOUSTIC 

C POWER. 

C PWR=THE ELECTRICAL POWER NECESSARY TO DRIVE THE SYSTEM 

C AT THE SPECIFIED INPUT CURRENT. 

C PUPWR= A POWER RATIO, ACOUSTIC TO ELECTRICAL. AN EFF I - 

C CIENCY . 

C 

C B . DATA DESCRIPTION. 

C FOR THE SELECTED INPUT DATA THAT IS HARD INPUT THE 

C FOLLOWING VALUES WERE EMPIRICALLY DETERMINED: M, L, Q, FO . 
C A, VOL, BL. S, SSUP, AND RE. RHO AND C WERE SELECTED FOR 
C HELIUM GAS AT 10 ATMOSPHERES AND A MEAN TEMPERATURE OF 300 
C DEGREES KELVIN. RM WAS CALCULATED FROM A THEORETICAL 
C EQUATION USING RHO , C S AND THE DRIVER RESONANT FREQUENCY. 
C RG WAS SET TO ZERO TO MODEL AN IDEAL POWER SOURCE. LE WAS 
C SET TO ZERO FOR EASE OF CALCULATION, BUT WAS TESTED DURING 
C RUNNING TO VERIFY THAT IT HAD LITTLE IF ANY EFFECT WHEN 
C INCLUDED. A CURRENT, I, OF ONE AMPERE WAS SELECTED FOR 
C EASE OF ANALYSIS. 

C SINCE YOU MAY WANT TO CHANGE MY INPUT DATA I WILL PUT 

C A COPY IN THE COMMENTS SO THAT YOU CAN RETURN TO IT IF 
C DESIRED. MY ANALYSIS HAS SHOWN THE BEST BACK VOLUME TO BE 
C 9.50E-8. AND THE BEST RM TO BE 1.535 KG/S. THIS VALUE FOR 
C RM WAS DETERMINED AFTER REVIEW OF PLOTS SHOWED HIGHER THAN 
C EXPECTED VOLTAGES AND IT WAS DETERMINED BY MEASUREMENT IN 
C THE LABORATORY. 

C 

c 

C DATA RM , M, SSUP , RHO, C/3 . 53E-3 , 0 . 0045 , 2 . 6E4 , 1 . 68 , 1000/ 

C DATA S, BL.Q.FO. A/1. 521E-3.8, 15,800,7. 13E-5/ 

C DATA VOL . RG ,LE,RE,I,L/2.715E-6,0,0,8,1,0.625/ 
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REAL LCOIL.K, IMZM.MAGU. IMZT ,M, L , LE , I , MAGV 
COMPLEX ZM , ZE , ZCOIL , ZTOT , V , U , ZTUEE , P 
P I = 3 . 141592654 

DATA RM , M , 3SUP , RHO .C/3 . 53E-3 , 0 . 0045 , 2 . 6E4 ,1.68, 1000 / 
DATA S.BL,Q,FO. A/1 . 521E-3 , 8 . 15 , 800 , 7 . 13E-S/ 

DATA VOL , RG , LE , RE , I , L/2 . 715E-6 . 0 , 0 , 8 , 1 . 0 . 625/ 

SVOL= ( A*A*RHO*C *C ) /VOL 
WO=2 * PI *FO 
ALPHA=WO/ (2*Q*C) 



CALCULATION SECTION: 

THE NEXT SIX LINES CREATE AN OUTPUT FILE FOR THE DATA 
THAT IS NECESSARY TO PRODUCE THE DESIRED PLOTS, AND PLACES 
HEADINGS ON THE DATA AS APPROPRIATE. 



OPEN ( UN IT =69 , FILE= ' FRIDGE OUTPUT A1 ’ , IOSTAT= IOS ) 

WRITE ( 69 , * ) ' MASS VOLUME' 

WRITE ( 69 , 1 ) M , VOL 

WRITE ( 69 , * ) ' FREQUENCY PIST VEL PRESS*PIST VEL PWR RATIO 

VOLTAGE’ 

1 FORMAT (E13.4.E13.4) 



THE FOLLOWING DO LOOP PROVIDES FOR THE CALCULATION OF 
THREE HUNDRED DATA POINTS FOR EACH OF THE DESIRED PLOTS, 
USING A FREQUENCY RANGE OF ONE QUARTER TUBE RESONANCE TO 
THREE HALVES TUBE RESONANCE. 



DO 10 J=1 , 300 

F=(F0/4)+( ( J-l ) * (FO/240) ) 



THE NEXT THREE LINES CALCULATE VARIOUS QUANTITIES, 
EXPLAINED ABOVE. THAT ARE FUNCTIONS OF FREQUENCY, FOR USE 
IN FUTURE CALCULATIONS. 



W=2 * PI * F 
K = W/C 

LCOIL=W* LE 



C THE NEXT NINE LINES CALCULATE THE REAL AND IMAGINARY 

C PARTS OF THE TUBE/DRIVER MECHANICAL IMPEDANCE AND COMBINE 
C THEM INTO A COMPLEX QUANTITY, ZM . 



REZM = RM+ ( ( ( RHO *C * S* ALPHA *COS(K*L) r SIN(K*L) ) + ( K * RHO * C 4 S 
*SINH ( ALPHA *L ) * COSH ( ALPHA * L ) ) ) / ( K* ( 1+ ( ALPHA/K) **2) * < ( ( ( 
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'* SIN ( K * L ) ) * *2) * ( ( COSH ( ALPHA *L) ) * *2) ) + ( ( (COS(K*L) ) '**2) * ( ( 

* SINH( ALPHA *L )) * *2) ))) ) 

IMZM= ( W*M ) - ( SSUP/W) - ( SVOL/W ) + ( ( ( RHO*C * S* ALPHA * SINH ( ALPHA 1 
1 *COSH ( ALPHA *L ) ) - ( K * RHO *C*S*COS ( K* L ) *SIN(K*L) ) )/(K* ( l+< ALP 

* /K) * *2) * ( ( ( ( SIN ( K*L ) ) ** 2 ) * ( ( COSH ( ALPHA * L ) ) **2) ) + ( ( (COS (K* 

* * * 2 ) * ( ( SINH ( ALPHA * L ) ) **2) ) ) ) ) 

ZM=CMPLX ( REZM , IMZM) 



THE NEXT THREE LINES CONVERT ZM TO AN ELECTRICAL 
QUANTITY AND THEN COMBINE ALL OF THE ELECTRICAL IMPEDANCES 
TO GET A TOTAL IMPEDANCE. 



ZE= ( BL* * 2 ) /ZM 
ZCOIL=CMPLX ( RE , LCOIL) 
ZTOT=RG+ZCOIL+ZE 



THE NEXT SIXTEEN LINES CALCULATE THE VALUES TO BE 
PLOTTED. AND PUT THEM IN THE OUTPUT FILE. AND THEN CLOSE 
THE FILE. 



V=I*ZTOT 
FORCE=BL* I 
U=FORCE/ZM 
M AGU =C ABS ( U ) 

REZT=REZM-RM 

IMZT=IMZM-(W*M) + (SSUP/W) + (SVOL/W) 
ZTUBE=CMPLX ( REZT , IMZT ) 

P= ( ZTUBE*U ) /S 

PU = 0 . 5 * ( REAL ( P*CON JG ( U ) )*S) 

PWR=0 . 5* ( REAL ( I *CONJG( V) ) ) 

PUPWR=PU/PWR 

MAGV=CABS ( V ) 

WRITE ( 69 , 100 ) F , MAGU , PU , PUPWR . MAGV 
100 FORMAT (E13.5,E13.5,E13.5,E13.5,E13.5) 
10 CONTINUE 
CLOSE (69) 



C PLOTTING SECTION: 

C THE NEXT TWO LINES INTERFACE WITH MY PLOTTING PROGRAM. 
C "GRAFRI FROTRAN Al", AND THE DISSPLA ROUTINE TO PLOT THE 
C DATA. 



PAUSE 'PRESS ENTER TO START PLOTTING ROUTINE’ 
CALL EXCMS ( ' EXEC DISSPLA GRAFRI’) 

STOP 

END 
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APPENDIX B. COMPUTER MODEL GRAPH PROGRAM 



PROGRAM GRAFRI 

C THIS IS LT. MICHAEL P. SUSALLA ’ S PLOTTING PROGRAM FOR 
C USE IN CONJUNCTION WITH "FRIDGE FORTRAN Al" IN "THE THER- 
C MOACOUSTIC HEAT ENGINE SIMULALT ION" . THIS PROGRAM MUST BE 
C COMPILED IN FORTVS PRIOR TO RUNNING "FRIDGE FORTRAN Al” IN 
C WF77 . THIS PROGRAM IS FOR USE IN DISSPLA AND CALLS VARIOUS 
C DISSPLA SUBROUTINES. THESE VARIOUS SUBROUTINES ARE USED 
C MORE THAN ONCE IN THE PROGRAM BUT WILL ONLY BE DESCRIBED 
C ONCE EACH. IN ORDER TO USE DISSPLA YOU MUST HAVE ARRAYS 
C OF DATA POINTS TO BE PLOTTED. MY ARRAYS ARE ALL DIMEN- 
C SIONED TO 300 POINTS SINCE THAT IS THE NUMBER CALCULATED 
C BY "FRIDGE FORTRAN Al". YOU SHOULD BE AWARE THAT IF YOU 
C CHANGE THE NUMBER OF POINTS TO BE CALCULATED, AND THE 
C NUMBER IS GREATER THAN 300, YOU MUST REDIMENSION THE 
C ARRAYS. TWO CHARACTER STRINGS ARE USED IN THIS PROGRAM. 
C NEW1$ IS USED IN RENAMING FILES TO THE FORMAT NECESSARY 
C FOR DISSPLA, AND DUMMY$ IS USED TO READ THE HEADINGS OUT 
C OF THE DATA FILE. ALL OF THE DATA TO BE PLOTTED IS IN THE 
C FORM OF REAL NUMBERS, SO ALL OF THE VARIABLES ARE MADE 
C REAL. D IS USED FOR THE FREQUENCY DATA, G IS USED FOR THE 
C PISTON VELOCITY DATA, H IS USED FOR THE ACOUSTIC POWER 
C DATA, O IS USED FOR THE ELECTROACOUSTIC EFFICIENCY DATA 
C AND, T IS USED FOR THE VOLTAGE DATA. THE XMAX AND YMAX 
C VARIABLES ARE INITIALIZED TO ZERO AND ARE USED TO CALCU- 
C LATE THE MAX VALUE OF EACH OF THE PARAMETERS AS THEY ARE 
C READ INTO THE ARRAYS. THIS COMPLETES THE INTRODUCTION AND 
C EXPLAINS THE FIRST FIVE LINES OF THE PROGRAM. 



DIMENSION D ( 300 ) ,G(300) ,H(300) ,0(300) ,T(300) 
CHARACTER *80 NEW1$,DUMMY$ 

REAL D , G , H , O , T , M 

DATA XMAX , YMAX1 , YMAX2 , YMAX3 , YMAX4/0 . , 0 . . 0 . , 0 . , 0 . / 



C THE NEXT SEVEN LINES OF THE PROGRAM RENAME MY DATA 
C FILE TO THE FORMAT NECESSARY FOR USE IN DISSPLA, OPEN THE 
C FILE, READ THE MASS AND VOLUME OUT OF THE FILE AND, READ 
C THE HEADINGS OUT OF THE FILE. 



NEW1$= ’ RENAME ’//' FRIDGE OUTPUT Al’//' FILE XYZ1 Al ’ 
CALL EXCMS ( NEW1 $ ) 

OPEN ( 69 , FI LE= ' XYZ1 ' ) 

READ ( 69 , ’ ( A80 ) ’ ) DUMMY$ 
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READ (69, 1)M, VOL 
READ ( 69 , ' ( A80 ) ' ) DUMMYS 
1 FORMAT ( E13 . 4 , E13 . 4 ) 



C THE FOLLOWING DO LOOP READS THE PLOTTING DATA INTO THE 
C ABOVE SPECIFIED ARRAYS FOR USE IN THE DISSPLA SUBROUTINES, 

C AND ALSO CALCULATES THE MAXIMUM VALUE OF EACH OF THE 
C ARRAYS FOR SETTING UP THE SCALES OF THE PLOT AXES. THEN 
C THE DATA FILE IS CLOSED. THESE VALUES ARE LEFT AS REAL 
C NUMBERS VICE INTEGERS BECAUSE THE DATA FROM THE PROGRAM 
C CAN VARY SIGNIFICANTLY WITH DIFFERENT INPUT VALUES, AND IF 
C THE MAX VALUE IS LESS THAN 10 AN INTEGER REQUIREMENT GIVES 
C A PLOTTING ERROR. 



DO 20 J = 1 ,300 

READ (69, IS, END=21)D( J) ,G(J),H(J),0( T ',T(J) 
15 FORMAT (E13.5,E13.5,E13.5,E13.5,E13.5) 

I F ( D ( J ) . GT . XMAX ) THEN 
XMAX=D ( J ) 

ENDIF 

I F ( G ( J ) . GT . YMAX 1 ) THEN 
YMAX1=G ( J ) 

ENDIF 

IF ( H ( J ) . GT . YMAX2 ) THEN 
YMAX2=H( J) 

ENDIF 

IF ( O ( J ) . GT . YMAX3 ) THEN 
YMAX3=0 ( J ) 

ENDIF 

IF ( T ( J ) . GT . YMAX4 ) THEN 
YMAX4=T ( J ) 

ENDIF 

20 CONTINUE 

21 CLOSE (69) 



C THE NEXT THREE LINES TAKE CARE OF ANOTHER ODDITY OF 
C THE DATA. IN THE CASE OF THE VOLTAGE PLOT. DEPENDING ON 
C THE INPUT DATA THE VOLTAGE WILL LOCK IN AT 8 VOLTS, AS 
C EXPECTED, BUT THIS IS ALSO THE MAX VALUE, SO YOU GET A 
C STRAIGHT LINE AT THE TOP OF THE PLOT MARGIN. THE NEXT 
C THREE LINES INCREASE THE MAX VALUE TO 9 VOLTS TO MAKE THIS 
C PLOTTING CASE MORE CLEAR. 



IF ( YMAX 4 . LE . 9 . 0 ) THEN 

YMAX4=9 . 0 

ENDIF 
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Q 0 0 0 0 



THE NEXT EIGHT LINES CALCULATE THE INCREMENTAL STEPS 
FOR EACH AXIS OF THE VARIOUS PLOTS. THE X AXIS IS THE 
SAME FOR ALL, AND IS LEFT IN INTEGER FORMAT FOR CLARITY. 
EACH AXIS IS SET UP FOR TEN INCREMENTS. THE DATA FILE IS 
THEN RENAMED TO THE DESIRED FORMAT FOR MY STORAGE. 



INC=INT (XMAX/10 . ) 
YINC1=YMAX1/10 . 
YINC2=YMAX2/1 0 . 
YINC3=YMAX3/1 0 . 
YINC4=YMAX4/1 0 . 
MAX= INT ( XMAX ) 
NEW1$= ' RENAME '//’ 
CALL EXCMS ( NEW1$ ) 



FILE XYZ1 A1 ' // ' 



FRIDGE OUTPUT A1 ' 



C THIS IS THE FIRST PLOTTING SECTION. IT DOES THE 
C INITIAL SET UP OF ALL OF THE CALLS THAT CAN BE COMMON AND 
C PLOTS THE VOLTAGE VERSUS FREQUENCY PLOT. CALL TEK618 SETS 
C UP THE REQUIRED MACHINE FOR THE PLOTTING AND LINKS YOU TO 
C THE PLOTTER. CALL BLOWUP IS USED TO MAGNIFY THE PLOT FOR 
C A BETTER PICTURE. CALL YAXANG PUTS THE AXIS NUMBERS 
C PERPINDICULAR TO THE PLOT AXIS TO PREVENT RUNNING INTO 
C EACH OTHER. CALL PAGE SETS UP THE PAGE SIZE IN INCHES. 
C CALL AREA2D SETS UP A TWO DIMENSIONAL PLOT AREA OF THE 
C SPECIFIED SIZE IN INCHES. CALL SCMPLX SELECTS THE LET- 
C TERING TYPE FOR THE HEADING. CALL MESSAG AND CALL REALNO 
C ALLOW ME TO PLACE A MESSAGE IN THE PLOT AREA AT A SPECIFIC 
C POINT SPECIFIED IN INCHES IN THE CALL. CALL XNAME AMD 
C YNAME ALLOW ME TO LABEL THE AXES. CALL HEADIN ALLOWS ME 
C TO PLACE A TITLE ON THE PLOT. CALL THKFRM AND FRAME ALLOW 
C ME TO PUT FRAMES AROUND THE PLOT AREA AND THE WHOLE PLOT 
C OF A SET THICKNESS. CALL GRAF SETS UP THE PHYSICAL ORIGINS 
C AND MAX VALUES AND THE INCREMENTS FOR EACH AXIS. CALL 
C CURVE TELLS THE MACHINE TO PLOT THE DESIRED ARRAYS OF DATA 
C POINTS. CALL ENDPL ENDS THAT SPECIFIC PLOT AND GIVES A 
C NEW PAGE FOR THE NEXT PLOT. 



CALL TEK618 
CALL BLOWUP (1.5) 

CALL YAXANG ( 0 ) 

CALL PAGE (11, 8. 5) 

CALL AREA2D ( 8 . ,6. ) 

CALL SCMPLX 

CALL MESSAG ( ’MASS=$' ,100,0.1,5.5) 

CALL REALNO (M, -2, ’ABUT’ , 'ABUT ' ) 

CALL MESSAG ( ' KG$ ' , 100 , ’ ABUT ’ , ’ ABUT ' ) 
CALL MESSAG ( ’ VOLUME=$ ' ,100,0.1,5.) 

CALL REALNO ( VOL , -2 , ' ABUT ’ , ’ ABUT ' ) 

CALL MESSAG ( 'M**3$' , 100, 'ABUT' , 'ABUT' ) 
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CALL XNAME( ' FREQUENCY ( IN HERTZ )$’, 100 ) 

CALL YNAME( ’VOLTAGE (IN VOLTS )$’. 1 00 ) 

CALL HEADIN( ’ VOLTAGE VS. FREQUENCY$ ’ ,100,2. , 1 ) 
CALL THKFRM ( . 02 ) 

CALL FRAME 

CALL GRAF (0.0, INC , MAX ,7.5, YINC4 , YMAX4 ) 

CALL CURVE (D,T, 300,0) 

CALL ENDPL(O) 



THIS SECTION PLOTS PISTON VELOCITY VERSUS FREQUENCY. 



CALL AREA2D ( 8 . , 6 . ) 

CALL MESSAG ( ’MASS=$' ,100,0.1,5.5) 

CALL REALNO ( M , -2 , ’ ABUT 1 , ’ ABUT ' ) 

CALL MESSAG ( ’ KG$ ' , 100 , ’ ABUT ’ , ’ ABUT ' ) 

CALL MESSAG ( ' VOLUME=$ ’,100,0.1,5.) 

CALL REALNO ( VOL , -2 , ’ ABUT ’ , ' ABUT ’ ) 

CALL MESSAG ( ’M**3$’ ,100, ’ABUT’ , ’ABUT’ ) 

CALL XNAME( 'FREQUENCY (IN HERTZ) $', 100 ) 

CALL YNAME( ’ PISTON VELOCITY (IN METERS/SEC )$', 100 ) 

CALL HEADINC ’ PISTON VELOCITY VS. FREQUENCY'S ’ ,100,2. , 1 
CALL THKFRM (.02) 

CALL FRAME 

CALL GRAF ( 0 . 0 . INC , MAX , 0 . 0 , YINC1 . YMAX1 ) 

CALL CURVE (D,G, 300,0) 

CALL ENDPL(O) 



THIS SECTION PLOTS ACOUSTIC POWER VERSUS FREQUENCY. 



CALL AREA2D ( 8 . ,6. ) 

CALL MESSAG ( ’ MASS=$ ’,100,0.1,5.5) 

CALL REALNO (M, -2. ’ABUT’ , ’ABUT’ ) 

CALL MESSAG ( ’ KG$ ’ , 100, ’ABUT’ , ’ABUT’ ) 

CALL MESSAG ( ’ VOLUME=$ ’ .100,0.1,5.) 

CALL REALNO (VOL, -2, ’ABUT’ . ’ABUT’ ) 

CALL MESSAG ( ’M**3$’ ,100, ’ABUT’ , ’ABUT’ ) 

CALL XNAME( ’ FREQUENCY ( IN HERTZ )$’, 1 00 ) 

CALL YNAME( ’ACOUSTIC POWER (IN N-M/S ) $ ’ , 1 00 ) 

CALL HEADIN( ’ACOUSTIC POWER VS. FREQUENCYS ’ , 1 00 , 2 . , 1 ) 
CALL THKFRM (.02) 

CALL FRAME 

CALL GRAF ( 0 . 0 , INC , MAX ,0.0, YINC2 , YMAX2 ) 

CALL CURVE (D,H, 300,0) 

CALL ENDPL(O) 



THIS SECTION PLOTS POWER RATIO VERSUS FREQUENCY. 
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CALL AREA2D ( 8 . ,6. ) 

CALL MESSAG ( ’ MASS=$ ' ,100,0.1,5.5) 

CALL REALNO ( M , -2 , ' ABUT ’ , ' ABUT ’ ) 

CALL MESSAG ( ' KG$ ' , 1 00 , ' ABUT ’ , ’ ABUT ’ ) 

CALL MESSAG ( 'VOLUME=$' ,100,0.1,5. ) 

CALL REALNO ( VOL , -2 , ’ ABUT ’ , ’ ABUT ’ ) 

CALL MESSAG ( ’M**3$' , 100, ’ABUT’ , 'ABUT' ) 

CALL XNAME( 'FREQUENCY (IN HERTZ )$’, 1 00 ) 

CALL YNAME( ’ ACOUSTIC PWR/ ELECTRIC PWR$’,100) 

CALL HEADIN( ' POWER RATIO VS. FREQUENCY$ ’ ,100,2. , 1 ) 
CALL THKFRM ( . 02 ) 

CALL FRAME 

CALL GRAF ( 0 . 0 , INC , MAX ,0.0 , YINC3 , YMAX3 ) 

CALL CURVE (D, 0,300,0) 

CALL ENDPL(O) 



CALL DONEPL ENDS THE PLOTTING SEQUENCE AND SENDS 
BACK TO CMS. 



CALL DONEPL 

STOP 

END 



YOU 
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APPENDIX C. STAR DRIVER DIMENSIONAL DRAWINGS 




Figure C-l . Assembly Drawing of NIB Driver with 
all Housing Parts and Test Lid. 
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DRIVER HOUSING - RESONATOR VIEW 

[Dimensions in inches ) 



Figure C-2. Driver Housing — Resonator View. 
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Figure C-3. Driver Housing--Side View Section AA 
Showing Electrical Leads, Bolt Holes, 
Microphone Penetration, and FET 
Follower Housing. 
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end FET locoilon. 
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ELECTRONICS PLUG 

(Dimensions In inches) 



Figure C-4. Feed-Through for Electrical Leads 
Shown in Figure V-5. 
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Figure C-5. Driver Housing--Side View Section BB 
Showing Pressure Gage and Helium 
Fill Ports. 
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Fill ports . 
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clearance holes For 

10/32 bolts at 
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DC PRESSURE TRANSDUCER 

(Dimensions In Inches) 



Figure C-6. DC Pressure Gage Feed- 
Through Plug Details. 
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DRIVER HOUSING - LID VIEW 



Figure C-7. Driver Housing--Lid View Showing 

Capillary, Microphone, and Follower, 
and Bolt Circles for Resonator 
Attachment, Pusher Plate, Pressure 
Lid and GAS Can Lid. 
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(Note: Plate is 3/10’ wide and 

and slopes 10° at outer 
edge . ) 



PUSHER PLATE (LID VIEW) 



(Dimensions in inches ) 



Figure C-8. Aluminum Pusher Plate. 
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(Note: O-ring groove 0.018’ wide 

and . 036 ’ deep : 

lid is 1/2’ wide ) 



PRESSURE LID (LID V I EW ) 

(Dimensions in inches) 



Figure C-9. Pressure Lid with Lead 
O-ring Groove. 
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O-RING GROOVE FOR 
PARKER 2-174 O-RING 

DEPTH § ; V 1 DTH 



Figure 
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C-10. Pressure and Helium Leak Test Lid. 
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REDUCER CONE 
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Figure C-ll . 



Reducer Cone. 
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APPENDIX D. STAR DRIVER LIST OF MATERIALS 



The driver consists of a modified JBL model 2450J neo- 
dymium-iron-boron magnet with iron pole pieces on either 
side. The diaphragm is titanium and the voice coil is 
aluminum . 

The housing, pusher plate, pressure lid, and reducer 
cone will all be machined from 6061-T6 aluminum. All open- 
ings in the driver housing will be sealed with 1/16" or 
1/32” diameter lead o-rings . The electrical feed-throughs 
will be sealed with Emerson & Cuming's STYCAST <_rM> 2850FT 
epoxy. All wire leads are copper. The steel bolts connecting 
the resonator to the driver are size 1/4-20. All other steel 
bolts are size 10-32. 
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